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ABSTRACT 


This thesis takes five FORTRAN IV programs from ''Com- 
puter Programs for Computational Assistance in the Study 
of Linear Control Theory" by Melsa and Jones and translates 
them into a microcomputer BASIC language to run on an inex- 
pensive microcomputer system. Three of the five programs are 
state variable programs. They are BASMAT for basic matrix 
manipulation, RTRESP for rational time response, and GTRESP 
for graphical time response. Two are transfer function pro- 
grams, FRESP for frequency response and RTLOC for root locus. 
A user's guide and example are included for each. 

A final example is used to demonstrate the utility of 
mie two transfer function programs as an aid to direct 


digital design in the w'-plane. 
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I. INTRODUCTION 


The purpose of this thesis is two fold. First, it 
is to show that standard computer programs useful in the study 
of linear control theory may be adapted to run on an inex- 
pensive home/microcomputer. To demonstrate this, five 
programs were chosen from [Ref. 1], three state variable 
programs and two transfer function programs. The three state 
variable programs are BASMAT, a basic matrix manipulation 
program, RTRESP, a rational time response program, and GTRESP, 
a graphical time response program. The two transfer function 
programs are FRESP, a frequency response program, and RTLOC, 
a root locus program. These programs as they appear in [Ref. 
1] are written in basic FORTRAN IV language to be run on 
a main frame computer system utilizing standard graphics 
Subroutines. These programs were modified and rewritten in 
a microcomputer BASIC language which is an interpreted 
language. Generally these programs are limited to systems 
Saeloth order. It is felt that this limitation is more than 
acceptable for the purposes of this thesis. Also no major 
effort has been made to analyze or improve the efficiency of 
the numerical methods used. An effort of this type is 
advised if these programs are to be modified for higher 


order systems. 





The second part of this thesis investigates the methods 
and relationships involved in direct digital design in the 
w'-plane. The two transfer function programs adapted to run 
Oh a microcomputer from the first part of this thesis, FRESP 
and RTLOC, are used to aid in this investigation. 

Section II summarizes the common problems and consi- 
derations involved in the translation of the five programs 
to be run on a microcomputer system. Section III describes 
the three state variable programs, BASMAT, RTRESP, and 
GTRESP, and gives an example of their use and output. 
section IV is Similar to section III and describes the two 
transfer function programs, FRESP and RTLOC. 

section V deals with the eae: Subsection A gives 
some background on the w'-plane and subsection B provides a 
Simple example using the two transfer function programs to 
compare the s and w' planes. Subsection C develops templates 
of some constant parameters in the w'-plane. Subsection D 
ends the section with a more involved example. 


section VI provides some conclusions and recommendations. 





II. TRANSLATION CONSIDERATIONS 


The five programs in this thesis, BASMAT, RTRESP, GTRESP, 
FRESP, and RTLOC, were translated directly from programs of 
the same name in [Ref. 1]. These programs were written in 
the basic FORTRAN IV computer language, a compiled language 
designed to run on a main frame or mini computer system. [In 
the first part of this thesis these programs are translated 
into a microcomputer BASIC language adapted to run on the 
microcomputer system described in Appendix A. 

In general the translation from FORTRAN IV to the micro- 
computer BASIC used posed no serious problems. Even though 
the BASIC language used is by necessity an abridged version 
of the BASIC language found on larger more expensive systems, 
it was extensive enough to provide the necessary commands 
for these programs. 

The input portion of the programs were rewritten to be 
interactive for convenient keyboard entry of problem para- 
meters eliminating input formatting errors. The programs 
requiring more extenSive input were given the added capability 
of saving and retrieving a problem description from a disk 
file for the user's convenience. 

Unlike FORTRAN IV the BASIC used does not have complex 


math capabilities. Therefore, mathematical operations on 





eemplex Quantities were programed separately for the real and 
imaginary parts. 

Another important difference in the languages affecting 
mais lation iS that BASIC has no provision for local variables. 
All variables in the main program and all subroutines are 
global. This caused some bookKeeping problems in translating 
Subroutines to prevent undesired side effects. Another compli- 
@avion was that FORTRAN IV considers the first four characters 
of a variable name for identification of that unique variable. 
The BASIC used considers only the first two. This further 
complicated the bookkeeping of variables and resulted in 
variable names being assigned just because they were different 
from the rest and with no relation to the quantity represented. 

The output routines were written to conform as close as 
possible to the FORTRAN IV version. Since the microcomputer 
system used is limited to eighty columns on eight and a half 
inch wide paper, provisions were made within the program to 
automatically switch to a condensed character font when 
necessary to output greater than eighty characters per line. 

The programs requiring graphical output created some 
Meegue problems. First, no standard library subroutines for 
graphics as used in the FORTRAN IV programs were available 
on the microcomputer system used. As a result all graphics 
routines had to be developed for the system used. Also the 
appropriate variables were redimensioned giving consideration 
to the resolution of the graphics available on the 


11, 





microcomputer system to optimize the program somewhat. 
secondly, although it is possible to mix graphics and text 
Sieche microcomputer system, it is not done in a straight for- 
ward manner and utilizes extra memory. For memory considera- 
tions and convenience it was decided to provide sufficient 
information for interpretation of the graphical output 

below each graph. It is felt that this method is satisfac- 
tory and creates little inconvenience to the user. To 

enhance the interpretation of the graphical output much 

time and effort was devoted to developing plotting routines 

to display the data in relation to axes labeled with tic-marks 
and boundaries of integer vice fractional values. 

A major obstacle that was overcome was the identification 
and correction of a memory management problem unique to the 
microcomputer system used. This problem affected only those 
programs requiring graphical output. The essence of the 
problem was due to the size of the programs and the number 
of variables used, parts of the program and stored variables 
were being over written in memory by the graphics routines. 
This problem was finally solved by making appropriate changes 
in the memory management scheme. For a more detailed explana- 


tion of the solution see Appendix A. 


ae 





III. STATE VARIABLE PROGRAMS 


The three state variable programs discussed in this sec- 
tion are modified versions of the programs of the same names 
found in [Ref. 1]. These programs may be used as tools for 
the analysis and design of linear control systems expressed 


in the following state variable form: 


PaGee eet) et UGE) 


GE). = ie GL) x (t)] 
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The Basic Matrix program (BASMAT) described in part A 
of this section is used to compute the determinant, inverse, 
characteristic polynomial, and eigenvalues of the square 
fiers A. iit 1s also used to determine the resolvent matrix 
ferea) Micdmorate  thansit1on Matrix exp(At). 

The Rational Time Response program (RTRESP), described 
in part B, 1S used to determine a closed form expression 
for the time response of a system. The input function r(t) 
is required to have a rational time response and no repeated 
elgenvalues are allowed in the combination of the system 
and input. 

The Graphical Time Response program (RTRESP), described 
in part C, is used to produce a graphical display of the time 


response of a system to an arbitrary input. 
ee 





RTRESP and GTRESP can be used to study open loop systems 
by letting K equal zero and unforced systems by letting r(t) 


equal zero. 


ies 





A. BASMAT 
1. Basic Matrix Program (BASMAT ) 
we Inbroduct ion 
When done by hand, matrix manipulations can be 
quite tedious and the chances of an error being made are 
great. In the study of linear state variable analysis, a 
computer program to do these manipulations is almost 
essential. It can do the necessary manipulations much more 
quickly and accurately and allow the user to devote his 
time and energy to design and analysis. 
b. Description of Program 
BAS het Oconto (takes a matrix A, 

mmmeonmpuLres the Geterminant of A (det A), the inverse of 
A oe), the characteristic polynomial (det(sJ-A)), and 
Seeenvalues of A (A; as i CR So Memes Oly Olina Mid trax 

Dis) = (s]-a)7? 
and the state transition matrix 

Dit) = exp(At) 
The state transition matrix is expressed as matrix coeffi- 
cients times the natural modes exp( \t) with the complex 
eigenvalues expressed as damped sine and cosine terms. The 
resolvent matrix is written as 

Dis) = adj(sJ-A)/det(sI-A) 
wiemuhe Matrix numerator, adj(sI-A), is output as matrix 
coefficients of powers of s so that it takes the form 
N-1 


ats a need ot S 


Cees ey IS 9 N 
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The BASMAT program interactively accepts input 
of a matrix and calls subroutines to perform the appro- 
priate calculations. The subroutines used are CHREQ, CHREQA, 
PROOT, DET, SIMEQ, and STMST. These subroutines are listed 
below with a brief description. 

CHREQ. This subroutine is used to determine the 
Characteristic polynomial det(sI-A), and the resolvent 
matrix (eI-A) iGo tet twee Ne Leverrler  ateori thm 


ow 


is used to compute the resolvent matrix 
Ee) = (s]-4) 


in the form 


N 
wae N- 


i=1 + 
Dis) = - DCs) 


~— 


femewn(s) 1s the characteristic polynomial det(sI-A). 

The coefficients of the characteristic polynomial 
are determined by subroutine CHREQA. 

CHREQA. This subroutine is called by CHREQ to 
Getermine the characteristic polynomial, det(sI-A), for the 
Beatrix. The principal-minor method is used. 

PROOT. This subroutine uses a modified Bair- 
stow method for root extraction to determine the roots of a 
polynomial with real coefficients. It is used to calculate 
Pee wecazenvalues of the A matrix. 

DET. This subroutine computes the determinant 
of a matrix. A gauss elimination method to place the matrix 
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muemoper triangular form 1s used. It is called by CHREQA to 
calculate sub-determinants. 

SIMEQ. This subroutine is used in the FORTRAN 
version of BASMAT in [Ref. 1: pp. 7,8] to determine the 
miverse Of the A matrix. In the BASIC version a machine 
language Subroutine [Ref. 2] is substituted for speed 
and convenience. 

STMST. This subroutine is used to compute the 
state transition matrix 

Pit) = exp (At) 
femme Matrix A. It uses the Sylvester Expansion Theorem. 
Peeenvalues of the A matrix must be provided. This routine 
can not handle duplicate eigenvalues. 
c. Program Translation Problems 
See section II, Translation Considerations. 
2. BASMAT User's Guide 
This program occupies about 7k bytes of memory and 
does not utilize the graphics pages, therefore it is not 
necessary to relocate the disk operating system. 
Push "return" after each input. 
STEP 1: 
SCREEN PROMPT: 
BASIC MATRIX PROGRAM 


PROBLEM IDENTIFICATION —-> 


16 





REMARKS : 
Label the problem. Type in any appropriate combina- 


tion of letters, numbers, and/or symbols, excluding 


fonmas. This input 1S limited to 255 characters. 


SLEEP 2: 
SCREEN PROMPT: 
PLANT ORDER? -> 


REMARKS : 


Enter the number that corresponds to the order of 


the A matrix. Maximum order allowed is ten. 


STEP 3: 
SCREEN PROMPT: 
INPUT PLANT MATRIX. 
ACrow,column )= 
REMARKS: 
Input the A matrix by typing in each element as 


prompted. The program will ask for each element 


beginning with the first row and going from left to 


Pont . 
STEP 4: 
SCREEN PROMPT: 
HARDCOPY? (Y/N) -> 
REMARKS: 


After this step 1s completed the program will output 


Mee lah «ator reterence-—) This will be followed by 


the inverse of the A matrix, the determinant of the A 
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Higwix, ule Mamrix coetficients of the resolvent 


matrix numerator, the characteristic polynomial 


coefficients, the eigenvalues of the plant matrix, 


and finally the elements of the state transition 


matrix. See Figure 1 for a sample output. 


Oe 


BASMAT Example 


This example will use the system matrix 


Ones 
A=} 0 0 
2 -3 - 


Or @ 


to demonstrate the use of the BASMAT program. It will refer 


to the steps described in the previous section, BASMAT 


User's Guide. 


Step 1. Enter "EXAMPLE <CR>'"' where <CR> = return. 
meeo 2. Enter "3 <CR>" 
eamoeme racer "© <CR> 1 <CR> 0 <CR> 0 <CR> O <CR> 1 <CR> +2 
mee —3 <CR> -3 <CR>" 
Step 4. Enter "Y <CR>" 

The resulting output is shown in Figure 1. These 
results are interpreted as follows: 

The plant matrix is shown just as it was entered with 
mows Horizontal and columns vertical. This is the format 


Mmeea for all matrix output. 


The inverse of the plant matrix is shown followed by 


the scalar determinant of the A Mites ie 
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The matrix coefficients of the resolvent matrix 
numerator are given as powers of s. The characteristic 
polynomial is listed as coefficients of powers of s. 

Meem this output, the resolvent matrix 

@(s) = adj(sI-A)/det(sI-A) 
Gee written. The characteristic polynomial, det(sI-A), 
is 24+3s+3s-+ts>. The first element of the resolvent matrix 


numerator is 34+3s+s- making the first element of the 


resolvent matrix 


34+35+8~- 


So) 
iM 04+35+3s-ts° 


The real and imaginary parts of the eigenvalues of 
the plant matrix are listed. Finally, the elements of the 
state transition matrix are given. The first element of 


the state transition matrix can be written 


2 OzoL Oot 


Wht) = 0.333e 4% 0.6676 cos 0.866t+1.15e Sane Oe oOG 1 
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BASIC MATRIX PROGRAM 
freee? IDENTIFICATION-EXAMPLE 
Meee AAR RAKRAKRARAZARARAAKARARARRE 
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THE PLANT MATRIX 
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=. oO 
1 
-¥ 


= aa 
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so 
4 
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Peeve lERMINANT OF THE PLANT MATRIX 


' 
bo 
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THE MATRIX 


THE MATRIX“ 
3 
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fae MATRIX 
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THE MATRIX 
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Figure 1 (Cont. ) 
THE CHARACTERISTIC POLYNOMIAL IN ASCENDING FPOBIERS 
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ioe nabtonal Time Response Program (RTIRESP) 


a. kntroduetion 
Frequently it is desirable to know the response 
of a system as a function of time. A computer program can 
determine this quicker and more accurately than by hand. 
b. Description of Program 
RTRESP determines the time response in closed 


form of the closed loop system 


X(t )=Ax(t )+bu(t) 
7 als 
UG) =a ae Slee | 


y(t)=c'x(t) 


Mieco aly initial conditions x(0) and input r(t) for t20. 
The system must have a rational Laplace transform R(s) with 
a pole-zero excess of at least one. [Ref. 1: pp. 11,12] 

The input r(t) is treated by forming a m 8_order 
dynamic system whose initial condition response is equal to 
mer fOr a specific set of initial conditions. This system 
is combined with the original system and then the complete 
response in closed form is determined from the subroutine 
STMST. The order of the combined system must be ten or 
less. 

Various primary and utility subroutines are used 
in RTRESP. The primary subroutines are listed below with a 


brief description: 
Ze 





CHREQA. This subroutine determines the charac- 
teristic polynomial, det(sI-A), for the matrix A using the 
principal-minor method. 

PROOT. This subroutine uses a modified Bairstow 
method for root extraction to determine the roots of a poly- 
nomial with real coefficients. It is used to determine the 
roots of the numerator and denominator of R(s) when they are 
entered in polynomial form. PROOT is also used to determine 
the eigenvalues of the combined system matrix (i.e. the roots 
of the characteristic polynomial determined by CHREQA). 

SEMBL. This subroutine determines the coeffi- 
cients of a polynomial from the roots of the polynomial. 

It is used to determine the polynomial coefficients of the 
numerator and denominator of R(s) when they are entered 

in the factored form. This subroutine together with PROOT 
provide the feature that R(s) may be entered in either 
factored or polynomial form. It will appear in the output 
in both forms. 

STMST. This subroutine computes the state transi- 
tron matrix 

Dit) amex ( AG) 
for a matrix A. It uses the Sylvester Expansion Theorem. 
mepenvalues of the A matrix must be provided. This routine 
can not handle duplicate eigenvalues, therefore it is 


necessary that the combined system and input have no 
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repeated eigenvalues. STMST is used to determine the state 
transition matrix of the augmented system. 
c. Program Translation Problems 
see section II, Translation Considerations. 
Zeer REOP User's Guide 
Push "return" after each input. 
STEP 1: 
SCREEN PROMPT: 
RATIONAL TIME RESPONSE PROGRAM (RTRESP ) 
THIS PROGRAM DETERMINES THE TIME RESPONSE OF A 
CLOSED-LOOP SYSTEM DUE TO SPECIFIED INITIAL 
CONDITIONS AND INPUT. SYSTEM MUST HAVE A RA- 
TIONAL LAPLACE TRANSFORM WITH A POLE-ZERO EXCESS 
OF AT LEAST ONE. 
PROBLEM IDENTIFICATION -> 
REMARKS: 
Label the problem. Type in any appropriate combina- 
tion of letters, numbers, and/or symbols, excluding 
commas and colons. This input is limited to 255 
characters. This is also the name used to save and 
retrieve the disk file containing the problem 
description. (See steps 2 and 3) 
OTEP 2; 
SCREEN PROMPT: 


WAS THIS PROBLEM DESCRIPTION PREVIOUSLY 
SAVED? (Y/N) -> 


REMARKS : 
This step provides the option of recalling a pre- 
viously saved vroblem description from the disk. 


24 





mor ts Gyped, the program looks for the problem 
description saved under the label typed in STEP 1. 
The program then runs tO completion using the 
retrieved program description. Ifan'"N" is 
typed, the program continues with STEP 3. 

pubP 3: 


SCREEN PROMPT: 


DO YOU WANT TO SAVE THIS PROBLEM DESCRIPTION? 
a Ne 


REMARKS : 
A positive response will save the problem descrip- 
tion, that will be entered in the following steps, 
to the disk under the label entered in STEP l. 
This problem may be retrieved during a later session 
by entering the proper label in STEP 1 and typing 
a 1neteeEP 2, 
SEpP 4: 
SCREEN PROMPT: 
ORDER OF THE SYSTEM? -> 
REMARKS: 
The total order of the system and the input must 
not exceed ten. 
peeP 5: 
SCREEN PROMPT: 
INPUT SYSTEM (A) MATRIX 


A(1,1) = 
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REMARKS: 
Enter the elements of the A matrix as prompted. 
The format is A(Crow,column) = 
STEP 6: 
SCREEN PROMPT: 
THE A MATRIX 
(display of the A matrix) 
ANY CHANGES? (Y/N) -> 
REMARKS: 
If an "'N" is typed, the program proceeds to STEP 9. 
To correct the matrix type "Y'"' and the program will 
proceed to STEP 7. 
STEP 7: 
SCREEN PROMPT: 


TYPE ROW,COLUMN OF THE ELEMENT TO BE 
CORRECTED —> 


STEP 8: 
SCREEN PROMPT: 
Pa Crow column ) = 
REMARKS: 
Enter the correct value. The program will return 
FO) STEP 6. 
SBP 9: 
SCREEN PROMPT: 
INPUT THE CONTROL (B) VECTOR 


Bi) = 
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REMARKS: 
ae mr ew oO VeEClLOD 1s Cnterea, a correction Sequence 
Similar to steps 6 through 8 is encountered. 
SIEP 10: 
SCREEN PROMPT: 
INPUT THE OUTPUT (C) VECTOR 
C(1) = 
REMARKS: 
ee o mL NCeeCmvectOon TS entered, a correction sequence 
Similar to steps 6 through 8 is encountered. 
So EP 11: 
SCREEN PROMPT: 
INPUT THE FEEDBACK COEFFICIENTS 
FEEDBACK COEFFICIENT (1) = 
REMARKS: 
After the feedback coefficients are entered, a 
correction sequence similar to steps 6 through 8 
is encountered. 
BiEP 12: 
SCREEN PROMPT: 
INPUT THE GAIN -> 
REMARKS: 


Enter the controller gain. 
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Spr 13: 
SCREEN PROMPT: 
INPUT THE INITIAL CONDITIONS 
XO(1) = 
REMARKS: 
After the initial conditions are entered, a 
correction sequence Similar to steps 6 through 
8 is encountered. 
STEP 14: 
SCREEN PROMPT: 
ENTER THE INPUT GAIN -> 
REMARKS: 
Pacer the gain Of the input function, RCs). 
pieP 15: 
SCREEN PROMPT: 


'KEY' = P FOR POLYNOMIAL COEFFICIENT FORM 
"KEY" = F FOR POLYNOMIAL FACTORED FORM 


INPUT 'KEY', ORDER FOR NUMERATOR POLYNOMIAL ~-> 
REMARKS: 
Choose the preferred method of entering the numerator 
BOltynemtal of the input function, RCs). 
STEP 16A: 
SCREEN PROMPT (assumes coefficient form chosen): 
INPUT POLYNOMIAL COEFFICIENTS 
INPUT COEFFICIENTS IN ASCENDING PWRS OF S 


IPUL COBRFICEENT OF S°O => 
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REMARKS: 


Enter the coefficient along with its algebraic 
Slgn as prompted. The program assumes that the 
coefficient of the highest power of s is one. 
STEP 16B (assumes factored form chosen): 
SCREEN PROMPT: 
INPUT POLYNOMIAL FACTORS 


LOR ME Ae LOR, 2 
REAL PART —> 


IMAGINARY PART -> 
REMARKS: 
Enter the factor with its algebraic sign. Do not 
enter the associated root. (e.g. For "(s-1)" enter 
"-j",) After the real part of the factor is entered, 
the program asks for the imaginary part, allowing 
renernouc Of quadradie factors. If a non=zero 
imaginary part is entered, the program automatically 
enters its complex conjugate. 
STEP 17: 


eeneeN PROMPT: 


P FOR POLYNOMIAL COBRFFICIENT FORM 
F FOR POLYNOMIAL FACTORED FORM 


cy 
ice 


INPUT ‘KEY', ORDER FOR DENOMINATOR POLYNOMIAL ~-> 


REMARKS : 


The denominator is entered just as described for 


the numerator in steps 15 and 16. 
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After this step is completed, the program will 
Output all input for reference followed by the 
time response of x(t) and the system output, 
Pace ) 
3. RTRESP Example 
It is desired to know the time response of the 


closed loop system 


Oa 0 

MCE. = Cit) eet CG) 
ue — ji | o a. 

Gwe otuGe = 10 O25) x(t )} 
y(t) = [1 1] x(t) 


memo oe Input function is given by 
Riese — 12.0 “(ill/s) 
fia che initial conditions are 
ScrEye— © 
Referring to the steps described in the pervious sec- 
tion, RTRESP User's Guide, this problem would be entered as 
Zot Lows ; 
Step 1. Enter "EXAMPLE <CR>" where <CR> = carriage 
igeiuirn . 
SueCerceanawscep 3. Enter "“"N “<CR>." 
peep 4. NEnter “2 <CR>." 
Sten oO,sepmter 'O <CR> 2 <CR> -1 <CR> -1l <CR>." 
Seep 6. Enter "N <CR>," 
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Pecpeo meen cere © <CR= I <€R>." Enter appropriate 
response for correction sequence. 

Step 10. Enter "1 <CR> 1 <CR>." Enter appropriate 
response for correction sequence. 

pecomnhnh., —Hnter “LL <CR>).0 <CR>.'"' Enter appropriate 
response for correction sequence. 

peepelZ2. Enter “DS <CR>." 

prep bo. Enter "“O <CR> O <CR>."' Enter appropriate 
response for correction sequence. 

Beco la. Enter "2 <CR>." 

Step. fo. Enter “"P,O <CR>." 

mcep LGA. Enter "I <CR>." 

preepeiy, Enter “F,1 <CR> O <CRo." 


The program will now run to completion resulting in the out- 
put shown in Figure 2. 

The output begins with the program name and problem 
Maentification. The A, B, and C matrices are shown for 
reference followed by the feedback coefficients and the 
the initial conditions and input 


controller gain. Next, 


gain are listed. The numerator and denominator polynomials 
of the input function are given in both polynomial coeffi- 
cient form and factored form regardless of the method in 
which they were entered. The polynomial coefficient form is 
Paver aS a list of coefficients from left to right in 
ascending powers of s. The highest power of Ss is always 
one. The polynomial factored form appears as a list of the 
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real and imaginary parts of each root of the polynomial. 
Notice that the numerator is only given in coefficient 
form since it is of zero order and no roots exist. 

The time response of the state x(t) is given in the 
form of vector coefficients of the various natural modes of 
the system and the input function. From this output 


x(t) is seen to be 


MiG 7es- (o" cos 1,71t = 1.7067. >” sin 1.71t + 1.67 


tI 


x(t) 


Emccee.: (| O. Sabjneabe ville 


tl 


Xo (t) 


The time response of the output y(t) is given as scalor 
coefficients of the same natural modes as x(t) so it is seen 


that 


MMEG7e° °° cos 1.71t + 4.19e°°°'°" sin 1.71t + 1.67 
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RATIONAL TIME RESPONSE 
Seeleieelt TODENTIFICATION -> EXAMPLE 
SE Pe A AAR AR AR AAR AA AAR RAR AAA RRARAARARRRARARRARES 


THE A MATRIX 


THE B MATRIX 


3 1 
THE C MATRIX 


1 1 


me@eweack COEFFICIENTS 


av a a x af ss NO Nw 4 x ye? 4% BA Me AA AA %e wee UF NS dK NP NA UO NF % a a SPF AH AA Me FP 8 Mt UY SN « x % a? 8° °° YAS Nt 
RKKAKRK ERE RRS RK ER KER KK KKK REE LKR ER RRA KE KA KER AREA KES 


INITIAL COMDITIONS, 4C@) 


3 0 


it 
ho 


RGAIN 


PemerearUR POLYNOMIAL OF R45) - ASCENDING FWRS OF 5S 


1 


Figure 2 RTRESP Output 
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Figure 2 (Cont. ) 


“ e 


SaVeMiINaATOR POLYNOMIAL GF RCS) - ASCENDING PHRS OF 
J i 


DENOMINATOR ROOTS ARE 
REAL PART IMAGINARY PART 


iy ) 
x x x se 4 4 x x a7 %F 4% uF 4, “Ff % ~ * » x a x a Ss x a S/S NS ae % NFs x VP SF oA AP NS Re et ot NP LP AM NA AE A fnl® 
KEKRK RRR KK RR RK RAK ALR AER RAK AR RA RK RAK RRR KRRARK ERR KES REESE 


7% 


fei] RESPONSE OF THE STATE X¢T) 
feet COEFFICIENT OF EXP(-1.7S23 TACOS(C 1. 21S? 1sés7T 
me OO6S64/ ) 


eee ix COEFFICIENT OF EXP<C-1.73) THSING 1.7189 1369)T 


meer Yl 7 ose See yale Oo 


eee ee tUR COEFFICIENT OF EXAPCHOT 
Peo olOoO6/ Peo 2o7 2a 10 


MU Ee a KA RAK RARARRARRARRAARAAARARKARAARARARARAR 
Smee RESPONSE OF THE OUTPUT YCT> 
SeeerrieleNy OF EAPC -1. 7S TACOS (1.7137 18657T 


meee 7O556/ 
Beeler Pil ieNT GF EXPC=$2. 7S, TASI OM IL. Pl SP 1 sss T 


$#.13284143 
epee COEFFICIENT OF EXP(8)T 
1.5555664/ 


4 aseh sv % x Ss av » as x x » a7 % x SF SN SYN svt x NAF 4AP SPH SF AFA BA KF HAF KS ., xg NF RNP SF AEA EEN 4, .* APF AZ BP 4S 
KK KKK KK RHA KKK RRR KR KR KR RK REFER KR RE K ESSE KE EEK EK KEN ERS 


Figure 2 
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Sem clRESP 
ieecraphical Time Response Program (GTRESP ) 
ais Wohweiceyrehbkens steve 
Knowing the response of a system as a function of 
time can be very helpful in the study and analysis of that 
system. The graphical display of this response can give 
much insight into a system's response characteristics. 
Db.  Deseription of Program 
GTRESP [Ref. 1: pp. 22-28] determines and 
graphically displays the time response of the closed loop 


system 


x(t) = Ax(t) + bu(t) 

i 
Gem Kfar) — eect) ] 
Glee ce x(t ) 


aeembO any initial conditions x(0) and input r(t). 

The basic purpose of this program is Similar 
to that of the RTRESP program described earlier. The 
difference between the two programs is that GTRESP will 
determine the time response for arbitrary input functions 
which may not have rational Laplace transforms but RTRESP 
requires a rational Laplace transform. Also GTRESP produces 
a graphical display of the time response instead of a 
closed form expression as in the RTRESP program. 

The GTRESP program uses a fourth-order Runge-Kutta 
numerical integration algorithm to calculate the time response. 
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The main subroutines used are CALCU, RUNGE, TRESP, and YDOT. 
Also various utility and plotting subroutines are used. 
CALCU. This subroutine is called by TRESP to 


determine the reference input r(t) and the control input 
- 
ie -K [r(t) - k x(t)] 


The reference input r(t) must be defined by the user by 
inserting the appropriate BASIC coding into the subroutine 
between line numbers 5010 and 5500. The reference input is 
represented by the variable R and the control input by the 
variable U. 

RUNGE. This subroutine is called by TRESP and 
contains the actual fourth-order Runge-Kutta integration 
algorithm. It must be executed four times for each inte- 
gration step. 

TRESP. This subroutine is the driving subroutine 
which calls the subroutines CALCU, RUNGE, and YDOT along with 
the necessary plotting routines. It calculates the time 
response of the closed loop linear system described by the 
input parameters and plots the desired variables. 

YDOT. This subroutine is called by TRESP to 
compute the derivative 
x(t )=Ax(t)+bu(t) 

It is designed to handle linear systems but can easily be 
modified to handle nonlinear and time varying systems which 
would give GTRESP a nonlinear and time varying capability. 
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ise 


STEP 1: 


c. Program Translation Problems 
see section II, Translation Considerations. 
GTRESP User's Guide 


Push "return" after each input. 


SCREEN PROMPT: 


GRAPHICAL TIME RESPONSE PROGRAM (GTRESP ) 


THIS PROGRAM DETERMINES AND GRAPHICALLY 
DISPLAYS THE TIME RESPONSE OF A CLOSED- 
BOOP sioti DUE TO SPECIFIED INITIAL 
CONDITIONS AND INPUT. 


PROBLEM IDENTIFICATION ~? 


REMARKS: 


SLEP 2: 


Label the problem. Type in any appropriate 
combination of letters, numbers, and/or symbols, 
excluding commas and colons. This input is 
limited to 255 characters. This is also the name 
used to save and retrieve the disk file containing 


the problem description. (See steps 2 and 3) 


SCREEN PROMPT: 


WAS THIS PROBLEM DESCRIPTION PREVIOUSLY 
SAVED! (Y/N) => 


REMARKS: 


This step provides the option of recalling a pre- 
viously saved problem description from the disk. 

If "Y" is typed the program looks for the problem 
description saved under the label typed in STEP 1. 
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The program then runs to completion uSing the 
retrieved program description. Ifan"N" is 
Eyped the program continues with STEP 3. 
STEP 3: 
SCREEN PROMPT: 


DO YOU WANT TO SAVE THIS PROBLEM DESCRIPTION? 
CaN) ==> 


REMARKS: 
The problem description will be entered in the 
following steps. 
A positive response will save the problem des- 
cription to the disk under the label entered in 
STEP 1. This problem may be retrieved during a 
later session by entering the proper label in STEP 
meand typing "'Y"' ain STEP 2, 
STEP 4: 
SCREEN PROMPT: 
ORDER OF THE SYSTEM? -—> 
REMARKS : 
The total order of the system and the input must not 
execcea Len. 
SrEP 5S: 
SCREEN PROMPT: 
INPUT SYSTEM (A) MATRIX 


A(1,1) = 
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REMARKS: 
Mier the appropriate element of the A matrix. The 
format is A(Crow,column) = 
STEP 6: 
SCREEN PROMPT: 
THE A MATRIX 
(display of the A matrix) 
ANY CHANGES? (Y/N) -> 
REMARKS: 
If an "N" is typed, the program proceeds to STEP 9. 
To correct the matrix type "Y" and the program will 
proceed to STEP 7. 
SlEP 7: 
SCREEN PROMPT: 


TYPE ROW,COLUMN OF THE ELEMENT TO BE 
CORRECTED —-> 


STEP 8: 
SCREEN PROMPT: 
A(Crow,column) = 
REMARKS: 
Enter the correct value. The program will return to 
STEP 6. 
STEP 9: 
SCREEN PROMPT: 
INPUT THE CONTROL (B) VECTOR 
B(1) = 
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REMARKS: 
mapeweune B vector is entered a correction 
sequence similar to steps 6 through 8 is 
encountered. 
Step 10: 
SCREEN PROMPT: 
Fieve ine OUTPUT (C) VECTOR 
G@s) = 
REMARKS: 
marian C VECLor 15 Entered a correction sequence 
Similar to steps 6 through 8 is encountered. 
STEP 11: 
SCREEN PROMPT: 
INPUT THE FEEDBACK COEFFICIENTS 
FEEDBACK COEFFICIENT (1) = 
REMARKS: 
After the feedback coefficients are entered a correc- 
tion sequence Similar to steps 6 through 8 is 
encountered. 
BEEP 12: 
SCREEN PROMPT: 
INPUT THE CONTROLLER GAIN -> 
REMARKS: 


Enter the controller gain. 


40 





STEP 13: 
SCREEN PROMPT: 
INPUT THE INITIAL CONDITIONS 
XOC1) = 
REMARKS : 
After the initial conditions are entered a correc- 
tion sequence similar to steps 6 through 8 1s 
encountered. 
STEP 14: 
SCREEN PROMPT: 


ENTER THE FOLLOWING PARAMETERS : 
NOUEe. (lF=224) 7 (DE*ER) =< L00 


EVE Eeak TEME CLZ) => 
REMARKS : 
Enter the initial time of the time interval of 
interest. Due to program constraints the initial 
and final times and the time step and frequency of 
output must be chosen so that they satisfy the 
igedkab LOn. 
(TF-TZ) / (DT*FR) =< 100 
omeP 15: 
SCREEN PROMPT: 
FINAL TIME (TF) <-> 
REMARKS: 
Peter toe final time of the Time interval of 
interest. 
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peenP 16: 
SCREEN PROMPT: 
TIME STEP (DT) ~-> 
REMARKS: 
Enter the time increment for each step. 
STEP 17: 
SCREEN PROMPT: 
FREQUENCY OF OUTPUT (FR) -> 
REMARKS: 
Enter an integer n and the program will print out 
data on every nth time step iteration. Ensure 
that the relation expressed in the remarks of STEP 
14 is satisfied. 
SeP 18: 
SCREEN PROMPT: 
YOU MAY PLOT UP TO 8 VARIABLES VS TIME. 


VARIABLE NUMBER VARIABLE NUMBER 


aC) 1 Xo i 8 
X2(T) 2 SoICT) 0 
MaCL) 3 ALG) tO 
X4(T) 4 ECT Le 
ob), o UCT) 12 
X6(T) 6 el) i 
CD) 7 R(T) 14 


HOW MANY VARIABLES TO PLOT? MAX=8 ~> 
REMARKS: 
You may plot up to 8 variables vs time. These 


variables will be plotted on the same plot in the 
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output. Enter the number of variables you want to 
appear on the plot. 
SreP 19: 
SCREEN PROMPT: 
TYPE THE VARIABLE NUMBER <CR> -> 
REMARKS: 
Type the number associated with the variable of 
interest. After carriage return <CR> is typed, the 
program will continue to prompt for another variable 
number until the number of variables that the user 
indicated in STEP 18 has been entered. 
S2HP 20: 
SCREEN PROMPT: 
YOU CHOSE THE FOLLOWING VARIABLES: 
(list of variable numbers entered in STEP 19) 
DO YOU WANT TO MAKE ANY CHANGES? (Y/N) -> 
REMARKS: 
To make a change in the variables to be plotted 
type a "Y" and the program will return to STEP 
ioweeert an 'N' is typed, the program will run to 
completion. 
3. GTRESP Example 
It is desired to determine and plot the time response 
of the error 
SG amt y(t) 


the input r(t), and the state variable Xo (t ) for the system 


43 





O a 0 


X(t) = GD) aes u(t) 
-1 -1 i 

Meee tt eet OunO.o) xCly 

ee eee eC) 


The time interval of interest is 
O<t<4 
The iteration step size (DT) is chosen to be 0.01 and printed 
output is desired every ten steps so FREQ equals 10. The 
initial conditions are 
x(t )=0 


and the input function is 


ae = oe treat 


O otherwise. 


Ge) 

In order to define this input function the following 
BASIC coding is inserted between line numbers 5010 and 5500 
in the subroutine CALC. 

eoz0 IF TZ > 1 GOTO 5040 

S030 R= 3: GOTO 5500 

9040 R = 0 

9500 REM END OF ROUTINE DESCRIBING RCT) 

Referring to the steps described in section 2, 
GTRESP User's Guide, this problem would be entered as 
follows: 

Step 1. Enter "EXAMPLE <CR>'" where <CR> = carriage 
Bevurn. 

Step 2 and Step 3. Enter "N <CR>." 

Seep 4. Enter "2 <CR>." 
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Step 
Step 
Step 
response for 
Step 
response for 
Step 
response for 
Step 
Step 
response for 
Step 
Step 
Step 


step 


given in step 14 of GTRESP User's Guide is met. 


mieeaint is (TF-TZ)/(DT*FR) = 


In this case 


a Enter 


G6. Enter 


o-.  kn cer 


correction 


LO; Enter 


correction 


Lie enter 


correction 


i Seen ter 


ioe. Enter 


correction 


14. Enter 


£0. Enter 


16. Enter 


ys Enter 


(4-0) /(.01*10 )=40 


Om Che ac hee On =e CRS. 

Ph eho. 

JORG beh. ss bnter appropriate 
sequence. 


Pie Ok- 1 “Cne..'' Enger appropriate 
sequence. 

Pe Cnh  Je- Ch. Sener apprOopriate 
sequence. 

oe CRS 

(Oe CR (= Che. Enter appropriate 
sequence. 

GO CR. 

ae CR 

Spode oC he. a 

LO -Chie. wr rEnsunvestaat Ene constraint 
That con- 


100 


so the constraint is met. 


Step 


Step 


or 


Se 


Enter 


Enter 


m3 <CR>. tf 


Mie Ch loco ne se GecOR-. 2 This causes 


the error e(t), the input r(t), and the state variable Xo (t), 


respectively, 


statement. 


to be plotted as desired from the problem 


45 





The program will now run to completion resulting in 
the output shown in Figure 3. 

The output begins with the program name and problem 
identification. The A, B, and C matrices are shown for 
reference followed by the feedback coefficients and the 
Seurroller gain. Next the initial conditions and time 
parameters are listed. 

The second page of output lists in tabular form the 
value of time t and the corresponding values of the output 
y(t), the control u(t), and all of the state variables. 

The user has no control over the variables output in this 
Porm . 

The third page of output is the graph itself. Below 

the graph is enough information to properly interpret the 


plotted data. 
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GRAPHICAL TIME RESPONSE 
PROBLEM IDENTIFICATION -> EXAMPLE 


x 


ot. 
ot 
ot 
ot 
+, 
+. 
24s 
a 
ot 
>t: 
3+ 
ot 
ot 
+ 


oy, 

% . a PP AF AP % 
Pee oh oe he ee 
Pi we wt oN S% PNY AN KN 


SREKLEKE SELES SE KEES SEEKER ESE EE SES 
ot oN ote OTN eee 2e nase ote oe oe wT Te PM a aoe CNR PN PH Le we 


THE A MATRIX 
8 ! 
-1 -1 
THE B MATRIX 
y 


ieee MATRIX 


Beeeemcr, COEFFICIENTS 
1 se 


SAIN = 5 


“™ 

% 
“St 

~ 


INITIAL - Cits] Db [ITIigts “¢ 
3 4 


INITIAL TIME = 3 FINAL TIME = 4 


“ a ees PF 2 1 PF UP NP ae Se ae a Sa Sat “2 SP RP UF EP UP OP UP NP ee ON ae SF MP Me TP AP gt OF Ff TP MP oe My PM Met Mo * 
~~ ~_ + + salen allies allie allies tal ae ee a ee a Seb Pe oe eo > eee eam 
FR TN PN wie wm PK NS MD FN PN P% FN Sie oe oN Me oN SN ooh TS oD ote oe OD hh AN Me EN PN OR SN oD F% te FD PR SR TN FR TN 


STEP = .81 FREQUENCY OF DUTPUT = 1% 
44x 


“% * ~ 


Figure 3 GTRESP Output 
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meeure 3 (Cont. ) 


GRAPHICAL TIME RESPONSE 
Seer IDENTIFICATION -2 EXAEeMPLE 


iT a. ey ACT) eee 
3 3 20 4 3 
J 2. 19944778 iW? e2ear 1h 111864791 2 .89333383 
a2 3.34846675 143947762 ,a720L206 3.454344 
=o 3.9263394 18.4766595 082798787 424374261 
4 2.01337897 Vaaee tara! 122749807 4.5837989 
1S §.28499782 9.463332? 1.6889279 4.57468 1112 
0599999998 = 691852439 3.37543244 2. 13399243 4.371621764 
699999996 = $99 142589 2 22599887 2.0v91396 | 3. 77328044 
299999994 = 46812831 1.31321188 2.973439926 JeveduZcdd 
OIPITIFIL «=—--§.. 28825128 1.14451316 3426194266 3.91838362 
799999939 = 584783255 1.93493374 gedoel7a74 2.005839 1 
1.99999999 3.462837346 -18.1932886 3.696938?7? - 8369653122 
1.19799999 = 1..69968535 =13. 1373297 3.9964866 - 1.83738 124 
1.29999998 = .Z58794762 -8.71392849 3.08641664 -3.99762187 
1.39999998 -. 748993826 -9.43799354 2.7636982 Saarecagrad 
1.499999938 = -1.4543845 ~2.7954868 1 2.020797 | “492728541 
1.57799998 -1.96407293 ~ 799975606 2.15/99317 ~4 99197969 
1.69999998 = -2.85899263 7148711883 1.77328419 Suncdecoorad 
1679999997 = -2..89293483 1.71656348 1.49998864 ~3.49944248 
1.37999997 = -2..81688576 2 .39893976 1.47518785 =sde47 co! 
Meeeri7??e)©=©=—- 4 ,.80192291 2.69973162 19963825 See Salads? 
2.07999996 -1.64667353 2.7468897 5479 19246 eee leddeec, 
memez77770 © - 1.4 1988919 Peet leocay 04799815 =lenardtice 
mee2777796 «= 118722972 2.49587785 » 1972873584 - 1,398 108331 
2637999995 =~. 964837042 222996893 9724495893 ~1,93648655 
mey7?776 =~ 700194579 1,73954372 - 8168649882 mi ever at. 
2.99999995 = -.974947458 1.6314458 ~ 0774387299 ~ 4979148738 
Erer77777) «= ~. 41685978! 1.33487654 -. 117184954 ~ 279748746 
279999995 =~. 284449821 1.85867972 salavdzales ~, 149429754 
Beer e77779 «=~ . 176816196 818784237 - 197449518 ~ 8299586731 
2099999993 = - 8928628427 sd7o121916 ~ 149985564 85592492 | 
3.39999994 = =, 8276993776 41322497 1 - 137599611 » 197891223 
3.19999994 =, 8198515288 «25422934? - 124743259 19379473 
3627999994 = ,8918356439 siAarearao -, 189485217 198442861 
3.37999994 = 8718853287 68992849217 ~ 9938872894 1448924! 
3649999993 =. 881588273 ~ 8119845112 ~ 8768934685 » 198398742 
3.59999993 = .8847985445 - 85793434 - 9416168885 » 144497353 
3.69999993 882838823 - 8877343728 ~ 8477441999 198382217 
3279999993 = 8772893216 - 184315886 s.tsogcdd 172 112852641 
3.89999993 =. 8694859387 -. 118744482 - 9251881781 894574 1163 
3.99999992 8684867868 ~, 1994696615 - 8156831488 07789499275 


Figure 3 
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PROBLEM IDENTIFICATION ~- EXAMPLE 
%X GRAPHICAL TIME RESPONSE <x 


ABSCISSA -> TIME AXIS 
ORDINATE -—> RESPONSE MAGNITUDE 
TIC MARKS SHOW INTERVALS OF UNITY 
THE PLOT FRAME LIMITS ARE: 
ABSCISSA, 8 TO 4 
ORDINATE, -5 TO 5 


SYSTEM RESPONSE 
VARIABLE SYMBOL 
Ae Yt 


ERRORX 
INPUTT 


Figure 3 
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IV. TRANSFER FUNCTION PROGRAMS 


The two transfer function programs discussed in this 
section are modified versions of programs of the same names 
fomnd in’ {Ref. 1]. 

The Frequency Response program (FRESP), discussed in 
part A of this section, determines and plots the frequency 
response of a transfer function over a specified range of 
frequencies. The output may take the form of rectangular 
Bode plots or a polar Nyquist plot or both as desired. 

The Root Locus program (RTLOC), discussed in part B, 
Serrculates and plots the root locus of a transfer function 
for a specified range of gains. It is also possible to 
Seearee a small rectangular section of the root locus for 


more detail. 
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pee PE RESP 
1. Frequency Response Program (FRESP) 
ae INtCLrOaduct iON 
The response of a system as a function of fre- 
quency iS avery important characteristic of that system. 
Some common ways of graphically displaying the frequency 
response of a system include the use of amplitude and phase 
Bode plots and Nyquist diagrams. These are valuable tools 
for system analysis. 
b. Description of Program 
The FRESP program [Ref. 1: pp. 105-113] is 
used to determine the frequency response of a rational trans- 


fer function G(s) of the form 


G(s) = KaES 

where 
CG) Mele 25S ee ee 
Eto) eal age > S)ean ns pecem ee) 


The output may take the form of Bode plots or a Nyquist dia- 
mma Oor both in addition to tabular data. 

FRESP gives the user the option of supplying 
discrete frequency values or allowing the program to linearly 
or logarithmically interpolate frequency values between two 


limit values. The complex number G(jGJ) can be computed for 


out 





ememearequency value. G(jG)) can be written in rectangular 
form as 

SOO ee Cl) et ACC) ) 
where R(G@)) and X(GJ)) are real values. The magnitude and 


phase of G(j GJ) can then be written as 


a 
2 


G(jQ@)) = [R*(W) + X*(W)1 
ewe) = areban X(G@))/RC(G)) 

The FRESP program uses the subroutines PROOT, 
PVAL, and SEMBL. The subroutines MAXI, PHNOM, GRAPH, and 
SPLIT that appear in the FORTRAN version do not appear in the 
BASIC version. They are incorporated into the main program 
and into various plotting subroutines written for the specific 
microcomputer system used. 

PROOT. This subroutine uses a modified Bairstow 
method for root extraction to determine the roots of a 
polynomial with real coefficients. It is used to determine 
the roots of the numerator and denominator of G(s) when 
they are entered in polynomial form. 

SEMBL. This subroutine determines the coeffi- 
cients of a polynomial from the roots of the polynomial. 
It is used to determine the polynomial coefficients of the 
numerator and denominator of G(s) when they are entered in 
the factored form. This subroutine together with PROOT 
provides the feature that G(s) may be entered in either 
factored or coefficient form and it will appear in the output 


ma oOoth forms. 
D2 





PVAL. This subroutine is used to evaluate a 
polynomial A(s) with real coefficients for 
s = PR + jPI 
PVAL is executed twice for each frequency value used, once 
for the numerator of the transfer function and again for 
the denominator. 

c. Program Translation Problems 

In addition to the programing considerations dis- 
cussed in section II, an output anomaly was traced to an 
apparent oversight in the subroutine PVAL. It was discovered 
mae a Zero Order numerator over a second order denomina- 
tor with a free s such as 

Ss Coe 20) 
was plotted as 

1/(s+10) 
The reason for this was that PVAL treated a zero order 
perynomial and a first order polynomial both as a first order 
polynomial which in this case resulted in the free s in 
the denominator being cancelled. The reason for this can 
be seen by examining the portion of FORTRAN code from PVAL 
in [Ref. 1: p. 164] repeated below: 

P=CMPLX(A(NN+1) ,0. ) 

DO 100 J=1,NN 

100 P=P*S+A(NN+1-J ) 

VR=REAL(P) 

In this portion of code NN represents the order 
Sr the polynomial being evaluated and the intention is that 
the DO loop be executed a number of times equal to the order 
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of the polynomial. But because of the nature of the DO loop, 
if the polynomial is of zero order (i.e. NN equals zero) 
the loop will be executed once just as if the polynomial 
was of first order [Ref. 3]. One way to correct this 
problem is to modify the section of PVAL as shown here: 
P=CMPLX(CACNN+1),0. ) 
Pee chin) tO. 110; 96 
90 DO 100 J=1,NN 
100 P=P*S+A(NN+1-J ) 
110 VR=REAL(CP) 
The essence of the modification is to skip the DO loop if 
the polynomial has order zero. Translating this modifica- 
tion into the BASIC version of FRESP solved the problem. 
2. FRESP User's Guide 
Push "return" after each input. 
EoeP 1; 
SCREEN PROMPT: 
FREQUENCY RESPONSE PROGRAM (FRESP ) 
THIS PROGRAM OBTAINS AND PLOTS THE FREQUENCY 
OF A RATIONAL TRANSFER FUNCTION OVER A SPECIFIED 
RANGE OF FREQUENCIES. BOTH RECTANGULAR BODE 
PLOTS AS WELL AS A POLAR NYQUIST PLOT CAN BE 
OBTAINED. 
PROBLEM IDENTIFICATION ~> 
REMARKS: 
Label the problem. Type in any appropriate combina- 
tion of letters, numbers, and/or symbols, excluding 


@emmas and colons. This input is limited to 255 


characters. 


o4 





STEP 2: 
SCREEN PROMPT: 
INPUT THE GAIN ~-> 
REMARKS: 
Peter the transfer function gain. 
STEP 3: 
SCREEN PROMPT: 


P FOR POLYNOMIAL COEFFICIENT FORM 
F FOR POLYNOMIAL FACTORED FORM 


eB Y * 
BicHy! 


INPUT 'KEY', ORDER FOR NUMERATOR ,A(S) -> 
REMARKS : 
Choose the preferred method of entering the numera- 
tor polynomial of the transfer function, G(s). 
STEP 4A (assumes coefficient form chosen): 
SCREEN PROMPT: 
INPUT POLYNOMIAL COEFFICIENTS 
INPUT COEFFICIENTS IN ASCENDING PWRS OF §S 
Een CORT ICILENT OF S70 => 
REMARKS: 
Enter the coefficient along with its algebraic sign 
as prompted. The program assumes that the coeffi- 
cient of the highest power of s is one. 
STEP 4B (assumes factored form chosen): 
SCREEN PROMPT: 
INPUT POLYNOMIAL FACTORS 


FOR FACTOR 1 
REAL PART ~> 
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IMAGINARY PART ~> 
REMARKS: 
meer the factor with its algebraic sign. Do not 
enter the associated root. (e.g. For "(s-1)" 
enter "-1".) After the real part of the factor is 
entered, the program asks for the imaginary part 
merowing for input of quadradic factors. Ifa 
non-zero imaginary part is entered, the program 
automatically enters its complex conjugate. 
STEP 5: 


SCREEN PROMPT: 


P FOR POLYNOMIAL COEFFICIENT FORM 
hb LOR? POLYNOMIALS PACTORED FORM 


he 
"KEY! 


ieee S hu ORDER FOR DENOMINATOR, BCS) => 
REMARKS: 
The denominator is entered just as described for 
the numerator in steps 3 and 4. 
STEP 6: 
SCREEN PROMPT: 
MINIMUM FREQUENCY VALUE ~-> 
REMARKS: 
Enter the minimum frequency value in radians per 
second of the frequency range of interest. 
STEP 7: 
SCREEN PROMPT: 


MAXIMUM FREQUENCY VALUE -> 
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REMARKS: 


STEP 8: 


Enter the maximum frequency value in radians per 


second of the frequency range of interest. 


SCREEN PROMPT: 


NUMBER OF FREQUENCY VALUES TO BE USED ~> 


REMARKS: 


STEP 9: 


Enter an integer from 1 to 200. Since the horizon- 
tal resolution capability of the graphics on the 
microcomputer system used is 191 pixels, any number 
of values greater than 191 adds no more detail to the 
Bode plots. For the Nyquist diagram the greater the 
number of values the greater the detail due to the 
nature of the plot. If a list of discrete fre- 
quency values will be supplied (i.e. option 1 in 

step 9) enter the number of frequencies in that 


tenes t . 


SCREEN PROMPT: 


O = LOGARITHMIC INTERPOLATION 
1 = DISCRETE VALUES SUPPLIED 
2 = LINEAR INTERPOLATION 


CHOOSE ONE -> 


REMARKS : 


Make a choice by entering the associated number. This 
step allows the user to choose the method used by the 
program to obtain discrete frequency values used to 


of 





generate data points. LOGARITHMIC INTERPOLATION pro- 
duces data points that appear equally spaced on a 
logarithmic scale between the minimum and maximum 
frequency values entered in steps 6 and 7. With this 
choice a Bode and/or Nyquist plot may be obtained. 
LINEAR INTERPOLATION produces data points that appear 
equally spaces on a linear scale. Only a Nyquist 
plot may be obtained with this choice. DISCRETE VALUES 
SUPPLIED allows the user to supply a list of 

discrete frequency values of interest. No plots 

may be obtained, however, with this option. Below 

is a table summarizing the available graphical output 
for each choice. 


Cr0OTCr BODE NYQUIST 


O YES YES 
1 NO NO 
2 NO YES 


STEP 10 (assumes choice 1 in step 9): 
SCREEN PROMPT: 
TYPE IN THE DISCRETE FREQUENCY VALUES, 
(number entered in step 8) VALUES NEEDED 
FREQ (1)? 
REMARKS: 
Enter the list of frequency values as prompted. When 
the last value has been entered the program will run 


Gomcompletion outputting data similar to Figure 4 with 


the exception of the plots. 
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STEP 11 (assumes choice O in step 9): 
SCREEN PROMPT: 
BODE PLOT? (Y/N) => 
REMARKS: 
The program outputs tabular data regardless of the 
choice made here. 
STEP 12 (assumes choice O or 2 in step 9): 
SCREEN PROMPT: 
NYQUIST SPLOT? (Y/N) ~> 
REMARKS: 
The program outputs tabular data regardless of the 
choice made here. 
After this step, the program runs to completion. 
See sample output in Figure 4. 
See kESP Example 
It is desired to know the frequency response of the 
ieansier function 
SCO. o eas 


G(s) = 


4+ 6s + a6 =a 5? 


The frequency range of interest is G = 0.1 rad/sec to 100 
rad/sec. A Bode plot of the magnitude and phase is desired 
but a Nyquist plot is not. The plots should be generated 
from 100 values of frequency logarithmically spaced from 


Ser to 100. 


Se 





Referring to the steps described in the previous 
section, FRESP User's Guide, this problem would be entered 
as follows: 

STEP 1. Enter "EXAMPLE CR "' where <CR> = carriage 
fmecurn. 

pepe 2. Enter "8 <CR>." 

Sone 3. Enter "F,1 <CR>." 

mob e 4B. Enter ''.5 <CR> 0 <CR>." 

pur 5S. Enter "P,3 <CR> 4 <CR> 6 <CR> 3 <CR>." 

mee 6.) Enter '.1 <CR>." 

Sree 7. Enter "100 <CR>." 

See 8. Enter "100 <CR>." 

PEEP Oo. Enter "'O <CR>." 

See Lis Enter "Y¥ <CR>." 

seeP 12. Enter "N <CR>." 

The program will then continue to completion pro- 
ducing the output seen in Figure 4. 

The output begins with the program name and problem 
identification. The transfer function gain is given followed 
by the transfer function numerator and denominator each 
listed as coefficients in ascending powers of s and then as 
the real and imaginary parts of the roots. 

The second page of output is tabular data. It 
consists of the radian frequency and the transfer function's 
corresponding real and imaginary parts, magnitude, and phase 
in radians and in degrees. Although 100 frequency values 
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were generated and their associated data points plotted, 
tabular data appears only for every other data point. To 
reduce unnecessary output and for formatting purvoses the 


following scheme is used for tabular data. 


NUMBER OF FREQUENCIES TABULAR DATA 
REQUIRED FROM STEP 8 PRINTED FOR 
O to 50 every freq. 
There 100 every other freq. 
EO tow loo every 3rd freq. 
ee jefe 21018 every 4th freq. 


The next two pages of output are the Bode plots for 
amplitude and phase. The plots are headed with the problem 
identification and type of plot and below each plot is 
Sierra crtent information to allow proper interpretation of the 
plotted data. Note that the phase angles are normalized to 
always remain between -180 and +180 degrees. 

This is the end of the output generated from the 
example as it was input since a Nyquist plot was not desired. 
A Nyquist plot is included as the last page of Figure 4 for 


the sake of completeness. 
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BeeQueNCY RESPONSE 
feel IDENTIFICATION —- EXAMPLE 
BAR AKARA AA KAA RKKRA RRR ARR AKRRRRRRARARRRRKKARAE 
GAIN=8 
Seer OR COEFFICIENTS - IN ASCENDING POWERS OF & 
ee) 
1 
NUMERATOR ROOTS ARE 
REAL PART IMAGINARY PART 


ae 2 8 


MemerilNATOR COEFFICIENTS - IN ASCENDING POWERS OF & 


4 
} 
3} 
i 
DENOMINATOR ROOTS ARE 
REAL PART IMAGINARY PART 
sl —ly Pee aed | 
aot lL. Sense 3 | 
al % 


s ~ % *, » s SP % * “, s * SP uF 4% x x sep % x .,” % PSs os Ma SF? SF uP NSN SP APM SH NA KP KAP NS a.” Af SPS x 5 * se? % x ? 
KKKKRKKKRKA KR KR RRR KEK KKK KKK ERR ESKER KR KR KEE RK KEK ERKES 


Figure 4 FRESP Output 
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Figure 4 


RADIAN FREQ. 


187226722 
» 123284674 
141747416 
. 162975884 
» 187381743 
»219443469 
0247707638 
» 2848835988 
32745491? 
0375493582 
»432876 13 

497782359 
0972236769 
097933228 
795463392 
869749888 
1.8000088! 
1.149757 

1.32174116 
1.91991189 
1.7475284 1 
2.98723382 
2.31812972 
2.469688781 
3.85385554 
eet 117177 
4.0378173 

464138888 
9133669927 
§. 13598734 
7 85488239 
9.1119884 

732683358 
1..7226724 
12.3284676 
14.1747419 
16.297388 
18. 7981745 
21.5443472 
24.7797639 
28 488359 1 
32.7454921 
37 6493587 
43.2876135 
497782365 
J? 2236776 
3.793323 
#3.8463342 
96. 9749819 
188.880882 


(Cont. ) PROBLEM 


REAL PART 


1.8178582?7 
1,82247187 
1.82756848 
1.83885834 
1,83874983 
1.86665859 
1.88698827 
1,11285044 
1.145751? 

1, 18695955 
1.23774461 
1.299837? 
1.371888 15 
1.49524 1893 
1,33795782 
1.62457894 
1.69230769 
1.71476725 
1.84366872 
1.9183691 
769794858 
338989997 

~. 134336879 

~ 9423312989 

- 99737737 

-. 473804416 

~ 486613557 

- 32894 128 

- 237108135 

~ 201847626 

~. 194626892 

- 118278988 

~ 898161962 

- 8685723775 

- 8528722615 

- 83958843 17 

~ 829941485 

~ 8226335036 

- 8171782888 
9138854894 

9 84425637E-83 
7 458269 18E-43 
-3.637817 168-83 
4.2659 1987E-83 
-3.22764824E-83 
-2.441946898E-93 
~ 1.34746 146E-83 
= 1.37 /6009 1E-83 
~1.95734283E-83 
=2 .79879835E-84 


[HAGINARY PART 


8516283748 
89863428 13 
8663172682 
8744898294 
8832923195 
9920337756 

- 188388178 
197136197 
2111131384 

» 1182146475 

» 181419968 

. 9386 198318 
84226767 13 

- 0288852383 

-. 110868565 

— 260681164 

~ 2461538466 

~ 732787826 
-1.87488971 
-1.49816618 
-1,7451883 
-1.7966987 
“409599564 
-1.16877489 

- 683954287 
sge700/79 
942379268 
seeeuy? 240g 

» 192546166 
9921933898 
8578 1862 
8389 127542 

- 8293788175 

- 8165835577 

- 8198538543 
S7eli247loee-es 
-4,665196 19E-83 
-3.8622238E-83 
-2,91116276E-93 
-1.32142716E-43 
-8,58521676E-84 
-5./78986204E-84 
-3.75450885E-44 
-2.46911164E-284 
-1.62396323€-84 
-1.84818684E-24 
=/a2od7 7are 0 
=4oae2da/ (7b 80 
~3.340846697E-85 
-2.0085 1985E-85 


MAGNITUDE 


1.81835782 
1,824 1588? 
1.03178223 
1.84 152594 
1.95424531 
1,87861542 
1.89152785 
1.11799565 
1. 1311288 
1.192865? 1 
1.24189277 
1,23 15368? 
1.37165954 
1.45256 188 
1.549877 11 
1.4453479 
1.75411684 
1,36478888 
1.76392254 
2.82289944 
1.79459981 
1.83654745 
1.98176172 
1.24324255 
738678461 
» 713374931 
2031884177 
596282748 
pap Or eosde 
221178189 
mdagrisre 
. 124515475 
98936654855 
879549 1932 
8931914145 
848199462 
8583826496 
822889267 1 
08172795174 
8138724492 


98824953 1E-83 
70402117 25E-93 
9.65830488E-83 
4.2738587E-93 

Ge 2o/oboenda 
2 444296 1SE-83 
1,04879722E-43 
1.3984 1956E-83 
1.0577888 1£-93 
3.88 129964E-84 
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POENT LR ICATTGN —- EXAMPLE 


PHASE (RAD) 


8587 189222 
68972928479 
8643257541 
9716873248 
8798378887 
8868884 184 
8928257474 
9999768558 
89649 1838 

8925986297 
88179463658 
9519899349 
82031978662 


- 8 143787124 


~ 8765349743 
159958477 
- 266292852 
-, 403858573 
Seats 149535 
— 799389724 
-1,86534135 
= 1.962897 19 
- 1.65704 786 
-1.91866312 
sen lgarcia ys 
-2.5058656 1 
-2.44192781 
-2,.58802432 
-2.53862774 
-2.71163854 
"2.77201 189 
Geedcor go! 
-2.36722885 
-2.78420789 
-2.936096?7 
~2.96344 186 
-2.78792446 
-3.8074061 
-3.02584745 
-3.84833486 
-3.85359432 
-3.84918255 
soe0/d0707) 
-3.83377726 
-3.89 132896 
-3.99787757 
ay Ooo cee 
-3. 18853252 
-3.11284111 
-3. 11658756 


PHASE (CE) 


2.99973123 
3.292846 
3.88559056 
4, 1973832 
453146892 
473258435 
deererareo 
304998249 
300490362 | 
960859542 
4.484274633 
3.99 124729 
1.76584889 
Seecuarae? 
439512253 
See loccue? 
mi dwceduces 
23, 1389419 
“23, 1325894 
-45../779981 
-$1.93979965 
-73.9424464 
~94 9418367 
a0 57a laed 
= 122.256 142 
-132.378578 
-137.368575 
-146.183583 
se loleloecos 
lao g000 
-158.553287 
Slolevarce! 
-184.27986? 
- 156. 464599 
~ 168226889 
Glicy eu i27 co 
171, 14398? 
Sli 2ca 137 
ei fgeceeal 
~174.193418 
Seder ocls? 
alvdeclicde 
- 178. 198883 
-176.487485 
-177 119788 
-177.499374 
eitedoels ay 
-178, 1858358 
ml (aasde(ec 
=178.367393 





Figure 4 (Cont. ) 


PROBLEM IDENTIFICATION - EXAMPLE 
®¥X¥ BODE PLOT (AMPLITUDE) xx 


ABSCISSA -> COMMON LOG OF FREQUENCY 

ORDINATE -> COMMON LOG OF AMPLITUDE 

MINIMUM FREQUENCY SHOWN ON ABSCISSA -1 RADIANS/SEC 
MAXIMUM FREGUENCY SHOKIN ON ABSCISSA 108 RADIANS/SEC 
Meamemrgol LIMITS OF BODE PLOT ARE +-86 DECIBELS 


_> 
— 
— 
—= 


Figure 4 
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Figure 4 (Cont. ) 


PROBLEM IDENTIFICATION - EXAMPLE 
XX BOOE PLOT (PHASE) Xx 


ABSCISSA -> COMMON LQG OF FREQUENCY 
ORDINATE -> PHASE (DEGREES) 

TIC MARKS SHOW MULTIPLES OF 9@ DEGREES 
MINIMUM FREQUENCY SHORIN ON ABSCISSA = .1 RADIANS/SEC 
MAXIMUM FREQUENCY SHOLIN ON ABSCISSA = 148 RADIANS/SEC 
MAXIMUM PHASE ON ORDINATE SCALE 70 DEGREES 
MINIMUM PHASE ON ORDINATE SCALE -18@ DEGREES 


Figure 4 
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Figure 4 (Cont. ) 


PROBLEM IDENTIFICATION - EXAMPLE 
XX NYQUIST PLOT xXx 


ABSCISSA -> REAL PART OF GC JW) 
ORDINATE -> IMAGINARY PART OF G¢JWD 
TIC MARKS SHOW INTERVALS OF UNITY 
AXES CROSS AT ORIGIN 


Figure 4 
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fee RLLOC 
1. Root Locus Program (RTLOC) 
ape Int roduct 10n 
When analyzing a system it is helpful to know the 
location of the closed loop poles in the s-plane. These 
closed loop poles are the roots of the characteristic 
equation and determine the basic characteristics of the 
transient response of a closed loop system. A root locus 
plot is a plot of the roots of the characteristic equation, 
usually as a function of the gain of the transfer function, 
and therefore it is a valuable tool for system analysis. 
b. Description of Program 
The RTLOC program [{Ref. 1: pp. 114-121] calcu- 
lates and plots the roots of the equation 
ier kG(s) = 0 
as a function of K. G(s) is assumed to be a rational function 


of the form 


_ Ns) 
G(s) = am 
and the root locus becomes the locus of roots of Di. (Ss) as K 


varies where 


D,(s) = Dts) + KN(s) 


RTLOC uses an algebraic plus linear progression 
scheme to vary K to give reasonable spacing of the roots. 
D,(s) is obtained for each value of K and the subroutine PROOT 
is used to calculate its roots. 
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The scheme used to calculate values of K assumes 
that K takes on only positive values. If K is to range 
through negative value, the value of smaller magnitude (less 
negative) must be used as the minimum value. The routine 
starts at the maximum (less negative) value and becomes 
increasingly negative until the lower limit is reached. 

If both positive and negative values are desired then two 
separate runs must be made with only positive and only 
negative values. 

The RTLOC program has a feature that allows the 
user to specify a range of s and &) values around an area 
of the root locus plot that is of interest. That area is 
then enlarged and more closely spaced values of K are generated 
giving more detail to the plot. 

The major subroutines used are PROOT and SEMBL. 
The subroutine SPLIT used in the FORTRAN version is replaced 
pyeaaditional coding in the main program and by various 
plotting subroutines. 

PROOT. This Subroutine uses a modified Bairstow 
method for root extraction to determine the roots of a 
polynomial with real coefficients. It is used to determine 
the roots of the numerator and denominator of G(s) when it 
is entered in polynomial form. It is also used to determine 
the roots of Ds) for each value of K generated. 

SEMBL. This subroutine determines the coeffi- 
cients of a polynomial from the roots of the polynomial. It 
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is used to determine the polynomial coefficients of the 
numerator and denominator of G(s) when they are entered in 
the factored form. This subroutine together with PROOT 
provides the feature that G(s) may be entered in either 
factored or coefficient form. It appears in the output in 
both forms. 

c. Program Translation Problems 

In addition to the programing considerations 
discussed in section II, two problems were encountered. 

When a range of negative gains is entered the program cal- 
culates only one value of K due to the logic of line 229 in 
the FORTRAN code. This works correctly only for positive 
ranges of gain. (Note that K is represented by the variable 
G in the computer program.) In the BASIC version coding 

was added to test the sign of the range of gains and modify 
the logic to correctly handle the negative case. 

The second problem was encountered when using 
gains of very small magnitude necessary in using this 
program for w'-plane analysis. The schemes used to generate 
values of K for a reasonable spacing of the roots when plotted 
worked fine for the magnitudes usually encountered in s-plane 
analysis but was not flexible enough to handle gains of much 
smaller magnitudes (e.g. 0 to 1E-3). To add this necessary 
flexibility the following FORTRAN lines used to generate 


waines of gain G, 
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22¢ G=1.15*(G+SIGNG*0.05) 
228 G=1.04*(G+SIGNG*0.02) 


were modified by replacing the constant values 0.05 and 0.02 
with the variables Dl and D2, respectively, where 


D1=ABS (GMIN-GMAX ) /700 
D2=ABS (GMIN-GMAX ) /1500 


so that the gain increment was a function of the range of 
gains of interest. 
2. RTLOC User's Guide 
Push "return" after each input. 
STEP 1: 
SCREEN PROMPT: 
ROOT LOCUS PROGRAM (RTLOC) 


THIS PROGRAM PLOTS THE ROOT LOCUS OF A 
PiSsCRESED SYSrTEM. 


PROBLEM IDENTIFICATION —> 
REMARKS: 
Label the problem. Type in any appropriate combina- 
tion of letters, numbers, and/or symbols, excluding 
commas and colons. This input is limited to 255 
characters. 
SrnP 2; 
SCREEN PROMPT: 
INPUT THE RANGE OF GAINS (MIN,MAX) -> 
REMARKS: 
Enter the range of transfer function gains of in- 
terest. Note that the minimum gain is the gain of 


lowest absolute value. 
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STEP 3: 
SCREEN PROMPT: 


P FOR POLYNOMIAL COEFFICIENT FORM 
F FOR POLYNOMIAL FACTORED FORM 


prcay 7 
"KEY ! 


INPUT 'KEY', ORDER FOR NUMERATOR, N(S) 7> 
REMARKS : 
Choose the preferred method of entering the 
numerator polynomial of the transfer function 
GCs). 
STEP 4A (assumes coefficient form chosen): 
SCREEN PROMPT: 
INPUT POLYNOMIAL COEFFICIENTS 
INPUT COEFFICIENTS IN ASCENDING PWRS OF s OR 
70 
weer eCORFEECIENT OF SO —> 
REMARKS: 
Enter the coefficient along with its algebraic sign 
as prompted. The program assumes that the 
coefficient of the highest power of s is one. 
STEP 4B (assumes factored form chosen): 
SCREEN PROMPT: 
INPUT POLYNOMIAL FACTORS 


FOR FACTOR 1 
REAL PART -> 


IMAGINARY PART ~—> 


Cu 





REMARKS: 


Enter the factor with its algebraic sign. Do not 
enter the associated root. (e.g. For "(s-1)" enter 
"-i{",.) After the real part of the factor is 
entered the program asks for the imaginary part 
mrowane for input of quadradic factors. If a 
non-zero imaginary part is entered the program 
automatically enters its complex conjugate. 

mee O: 


SCREEN PROMPT: 


P FOR POLYNOMIAL COEFFICIENT FORM 
F FOR POLYNOMIAL FACTORED FORM 


"KEY! 
‘KEY ! 


INPUT 'KEY', ORDER FOR DENOMINATOR, D(S) ~? 
REMARKS: 
The denominator is entered just as described for the 
numerator in steps 3 and 4. 
SreP 6: 


SCREEN PROMPT: 


WOULD YOU LIKE TO LOOK AT ONLY A PART OF THE 
ROOT LOCUS? (Y/N) 7> 


REMARKS: 


Enter an "N'"' for the option of viewing the root locus 
for the entire range of gains entered in STEP 2. 
Bacer a “Y" for the option of viewing only a portion 
of the root locus defined by minimum and maximum 
values of sigma and W of interest to the user. 
Usually, it 1s best to view the entire root locus 
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first then decide if any section needs to be 
enlarged to reveal more detail and rerun RTLOC 
meeview Only that portion. 
STEP 7 (assumes "Y" entered in step 6): 
SCREEN PROMPT: 
ENTER SIGMA MIN, SIGMA MAX -> 
REMARKS: 
Enter the minimum and maximum values of Sigma (i.e. 
the real axis) of the portion of the root locus 
plot of interest. 
STEP 8 (assumes "Y" entered in step 6): 
SCREEN PROMPT: 
ENTER OMEGA MIN, OMEGA MAX ~-> 
REMARKS: 
Enter the minimum and maximum values of omega (i.e. 
the imaginary axis) of the portion of the root 
locus plot of interest. 
maenP 9; 
SCREEN PROMPT: 
PRINT OUT OF GAIN DATA? (Y/N) -> 
REMARKS : 
This gives the user the option to suppress the print 
out of the gain data if it is not needed. This saves 
time and paper if only the root locus plot is desired. 


After this step the program runs to completion. 
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oe RILOC Example 
It is desired to know the root locus for the open 


loop transfer function G(s) given by 


oa eS 


G(s) = 5 
s(8 + 9s + Ss") 


for gain K of from 0 to 30. Referring to the steps des- 
cribed in section 2, RTLOC User's Guide, this problem would 
be entered as follows: 

STEP 1. Enter "EXAMPLE <CR>" where <CR> = carriage 
return. 

Pepe 2. Hnter "0,30 <CR>." 

pee 3. Enter "“F,1 <CR>." 

eee 4B. Enter "1.2 <CR> O <CR>." 

Pitmeowerrncer 7,3 <CR> 0 <CR> 0 <CR> 1 <CR> O <CR> 
Seen © <CR>." 

STEP 6. Enter "N <CR>." 

eleP 9. Enter "Y <CR>." 

The program will then run to completion producing 
the output seen in Figure 5. 

The output begins with the program name and problem 
identification. Next the transfer function numerator and 
denominator are listed as coefficients in ascending powers 
of s and then as the real and imaginary parts of the roots 
(i.e. open loop zeros for the numerator and open loop poles 
for the denominator). The minimum and maximum values of gain 


are listed next. 
TA 





Peeene Option tO view only a part of the root locus 
is selected, this fact is stated next followed by the ranges 
of sigma and omega that are chosen. (See Figure 6) 

The next several pages of output contain the gain 
data which may be suppressed with the appropriate response in 
step 9. The data point number is listed with the value of 
gain and the real and imaginary parts of the corresponding 
roots of the open loop system. 

The last page of output contains the root locus plot 
itself. It begins with the problem identification and 
heading and ends with sufficient data listed to allow proper 
interpretation of the plot. 

Figure 6 shows the output from the same example 
problem with the gain data suppressed and the option to 


view only a part of the root locus chosen. 
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PROBLEM IDENTIFICATION - EXAMPLE 
XX ROOT LOCUS PLOT xXx 





ABSCISSA -> REAL (SIGMA) AXIS 
ORDINATE -> IMAGINARY (OMEGA) AXIS 
TIC MARKS SHOW INTERVALS OF 1 
THE PLOT FRAME LIMITS ARE: 
ABSCISSA, -8 TO 1 
ORDINATE, -4 TO 5 


Figure 5 
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ROOT LOCUS 
PROBLEM IDENTIFICATION - EXAMPLE WITH OPTICN 
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PROBLEM IDENTIFICATION - EXAMPLE WITH OPTION 
XX ROOT LOCUS PLOT xXx 





ABSCISSA -> REAL (SIGMA) AXIS 
ORDINATE -> IMAGINARY (OMEGA) AXIS 
TIC MARKS SHOW INTERVALS OF 1 
THE PLOT FRAME LIMITS ARE: 
ABSCISSA, -2 TO 1 
ORDINATE, -! TO 2 


Figure 6 


83 





V. w'!-PLANE ANALYSIS 


Digital control systems are becoming more and more common. 
Digital control laws have unique characteristics that can 
only be approximated by uSing classical techniques in the 
Sonu inuous s domain. 

In the w' domain all analog control system design tech- 
nology transfers completely for digital control system design. 
An important advantage of the w' domain is that non-minimum 
phase effects of the sampling and data-hold operations and 
of sampling rate can be directly accounted for without 
approximation while using conventional frequency domain 
design and analysis tools such as root locus and Bode 


plots. [Ref. 4] 


84 





A. BACKGROUND 
This s domain is used for continuous system analysis. 
When a digital system is considered the z, w, or w'domain 


must be used. These domains are related as shown below. 


ST 
Z=e 
zZz—-1 ST 
w= = tanh — 
ZtlL 2 
w' = (2/T) w 


where T is the sampling period and 


tanh x = 


In an s-plane root locus plot the region of stability is 
the left half plane, that is all roots with negative real 
Meros. This region of stability 1S mapped into a unit circle 
with its center at the origin in the z-plane. By use of 
the bilateral transformation shown above, the stability region 
of the z-plane is mapped back into the left half plane to 
form the w-plane. The w'-plane takes it a step farther by 
multiplying the w-plane by the factor 2/T where T is the 
Sampling period. This gives the w'-plane the property 
that not only is the region of stability the left half plane 
as in the continuous s-plane but w' approaches s as the sam- 
pling period T approaches zero [Ref. 5]. 
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B. TRANSFER FUNCTION APPROACH USING THE w'-PLANE 


In this section the second order system 






10 


s(st+10) 






is analyzed using the two transfer function programs dis- 
cussed in part IV of this thesis. It is then converted 
to the w'-plane and analyzed for periods of .001 seconds, 
.O1 seconds, and .1 seconds using the same two transfer func- 
tion programs. These results are used to gain insight into 
w'-plane analysis. 

The system is converted to the w'-plane by adding a 
sampler and a digital to analog converter in the form of a 


zero order hold. This modified system is shown below. 


i 10K 
SS Sete ae eee pT 
= S(s+10) 





The new open loop transfer function is transformed to 


the z-plane and then to the w'-plane as follows: 


ST 4 


meso 10K(l-e ~ ) 


s“(st10) 
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ma 


m.! 4 Tz Glacmae 7 
Gz) = 1OK(1-z ~) — —_— o ae ss TOT 
10 (z-1) 10(z-1)(z-e ) 
c ees 
= ae a 
z-1-10(z-e 197) 


For the period 
T=.001 


G(z) reduces to 


4.98337 x 107° K(z+.996681) 
ey ee (1) 
(z-1)(z-.990049834) 


Semrlarly, for 


T=.01 
4.837418 x 107* K(z+. 9672185) 
ee 
(z-1)(z-.904837418) 
id tor 
m=. 1 
03678794412 K(z+.718281827) 
oa 


(zZ-1)(z-. 3678794412 ) 


Now G(z) is converted to the w' domain using the rela- 


tionship 
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m+ T/2 w' 


me— T/2 w' 
For period 
T=.001 


1 + .0005 w' 


1 - .0005 w' 


Substituting into equation (1) gives 


6 1+.0005w' 
WegSoor x LO K ——_——————— +  , 996681 
1-.0005w' 
G(w') = 
1+.0005w' 1+.0005w' 
———— ——  - | —————— - .990049834 
1-.0005w' 1-.0005w' 


which reduces to 


4.156878 x 10 K (wt + 1202819) (w' - 2000) 


| 
Se een oy) 


Note that the gain is 1.0K which is the same as the s domain 


Se. oimilarly, for 


m=. O11 
1+.005w' 

) 
1-.005w' 
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-4.1625027 x ome K(w'+12002.00418 )(w!-200 ) 


G(w!) = —_— eee) 
w' (w! + 9.991674985) 


ama for 
T=, 1 
1+.05w' 
7; rere 
1-.05w' 


-.0037882843 K(w' + 121.9858708)(w'-20) 
a (45) 


AV en) vet ce Lats ses) 


Note again that the gain is 1.0K also in equations (3) and 
(4). 
1. Frequency Response 
a. sS-Plane 
The open loop transfer function of the system 
in the s domain is 
LOK 


G(s) a 
s(st+10) 


The frequency response of G(s) with a gain K of 5 is described 
by the output of the FRESP program shown in Figure 7. From 
the Bode plots it is seen that the gain margin is infinite 


and the phase margin is 65.1 degrees. 


Se 





b. w't-Plane with a Period of .001 Seconds 
The open loop transfer function for the equiva- 
lent sampled system is represented by equation (2) for a 
@ertoa fT of .0O01 seconds. This transfer function is 
entered into the FRESP program as if w' were an s. The 
resulting output in Figure 8 is interpreted as described 
earlier in the thesis keeping in mind that all references 
Boemseare actually referring to w'. From this output it 
is seen that the gain margin is no longer infinite as it 
wast the continuous case but it is still quite high at 50 
dB. The phase margin has dropped slightly from 65.1 degrees 
to 64.6 degrees. 
c. w'-Plane with a Period of .01 Seconds 
If the sampling rate of the system is decreased 
so that the sampling period is .01 seconds, then the open 
loop transfer function is represented by equation (3). 
Inputting the transfer function into FRESP results in the 
output shown in Figure 9. From the Bode plots the gain 
Marein 1s found to be 32 dB. This is down from 50 dB for 
the case of a period of .001 seconds indicating a decrease 
in stability. The phase margin is 64.5 degrees which is also 
lower than the previous case but only by one tenth of a 
degree. 
d. w'-Plane with a Period of .1 Seconds 
If the period iS again increased to .1 seconds, 
the open loop transfer function is represented by equation 
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foe Entering this transfer function into FRESP produces 
mires Oucput seen in Figure 10. From this output the gain 
margin is found to have decreased to 13.2 dB and the phase 
margin to 53.7 degrees. 
e. Summary of Frequency Response Results 
The table below is a brief summary of the gain 


and phase margins found in each case. 


CASE GAIN MARGIN PHASE MARGIN 
s-plane infinite So. deo. 
awe, 6 6T=.O001 Oued 64.6 deg. 
mw, 1=.01 32 GB 64.5 deg. 
we, T=.1 Powe. a5 93.7 deg. 


It can be seen from this table that the con- 
tinuous case is the most stable. The sampled cases become 
less stable as the sampling period increases. It is also 
noticed that the gain margin iS more sensitive to changes 
in the sampling period than the phase margin. 

fAeeeenOoot Locus 
Now a comparison is made using the root locus program 
uSing the same cases used above. 
a. s-Plane 
If the open loop transfer function 
10K 
G(s) =-——-— 
s(s+10) 
is entered into the RTLOC program, the result is the output 
in Figure 11. The root locus plot shows that the system 
never becomes unstable at any gain. 
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b. w'-Planewith a Period of .001 Seconds 
If the equivalent sampled transfer function with a 
sampling period of .001 seconds represented by equation (2) 
is entered into RTLOC, the output will be that found in 
Figure 12. From the root locus plot it can be seen that 
there is a slight tendency for the plot to curve toward the 
right half plane as it moves further from the real axis. The 
tendency is so slight that it will take a very large gain 
to drive the system unstable. 
c. w'-Plane with a Period of .01 Seconds 
If the sampling period is increased to .01 seconds 
the transfer function is represented by equation (3). 
emeering this transfer function into RTLOC results in the 
output seen in Figure 13. It can be seen that the tendency 
for the curve to bend toward the unstable right half plane 
1s increased indicating that a lowér value of gain than in 
the previous case will drive the system to instability. 
d. w'-Plane with a Period of .1 Seconds 
When the sampling period is further increased to .1 
seconds represented by equation (4), the resulting RTLOC 
output is that seen in Figure 14. In the root locus plot 
for this case the tendency to become unstable is much more 
pronounced. Here the plot bends into the right half plane 
within the limited boundaries of the portion plotted. It 


becomes unstable for values of gain K greater then 24.6. 
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e. Summary of Root Locus Comparison 


It can be seen that the continuous system 


10K 
= 
S(s+10) 
is always stable. When this system is sampled it can be 
seen by equations (2), (3), and (4) that a zero is added 
in the right half plane making it a non-minimum phase 
system. It can also be seen that by decreasing or increasing 
the sampling period the distance of this zero from the 
Origin increases or decreases respectively. The closer 
this zero is to the imaginary axis, the greater effect it 
has on bending the root locus into the right half plane. 


This effect can be seen by examining the root locus plots 


mame toOures 11 through 14. 
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Figure 7 (Cont. ) 
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XX BODE PLOT CAMPLITUDE) Xx 


ABSCISSA -> COMMON LOG OF FREQUENCY 

ORDINATE -> COMMON LOG OF AMPLITUDE 

MINIMUM FREQUENCY SHOWN CN ABSCISSA -1 RADIANS/SEC 
MAXIMUM FREQUENCY SHOMIN ON ABSCISSA 186 RADIANS/SEC 
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Figure 7 (Cont. ) 


PROBLEM IDENTIFICATION ~ EX1 S-PLANE 
XX BODE PLOT (PHASED XxX 


ABSCISSA -> COMMON LOG OF FREQUENCY 
ORDINATE -> PHASE (DEGREES) 

TIC MARKS SHOW MULTIPLES OF 96 DEGREES 
MINIMUM FREQUENCY SHOWN ON ABSCISSA = .1 RADIANS/SEC 
MAXIMUM FREQUENCY SHOWN ON ABSCISSA = 1468 RADIANS/SEC 
MAXIMUM PHASE ON ORDINATE SCALE = -986@ DEGREES 
MINIMUM PHASE ON ORDINATE SCALE = -186 DEGREES 
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Peeure 8 (Cont. ) 


PROBLEM IDENTIFICATION - EX1 W’-PLANE T=.661 EXTENDED 
5% BUDE PLOT (AMPLITUDE) kx 





ABSCISSA -> COMMON LOG OF FREQUENCY 

ORDINATE -> COMMON LOG OF AMPLITUDE 

MINIMUM FREQUENCY SHOMIN ON ABSCISSA = .1 RADIANS/SEC 
MAXIMUM FREQUENCY SHOWN ON ABSCISSA = 1086 RADIANS/SEC 
Pig@eeweneOe LIMITS GF BODE PLOT ARE +-166 DECIBELS 
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PROBLEM IDENTIFICATION - EX1 W’-PLANE T=.G61 EXTENDED 
Ao OWE See iG e hese jar 





Peeeisen —> COMMON LOG OF FREQUENCY 
GROINATE -> PHASE (DEGREES) 

toe MARKS SHOW MULTIPLES OF 96 DEGREES 
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MINIMUM PHASE ON GRDINATE SCALE piled DESREES 
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PHASE (RAD) 
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Meeure 9 (Cont. ) 


PROBLEM IDENTIFICATION - EX1 W’-PLANE T=.61 
| XX BODE PLOT CAMPLITUDE) xx 


ABSCISSA -> COMMON LOG OF FREQUENCY 

ORDINATE -> COMMON LOG OF AMPLITUDE 

MINIMUM FREQUENCY SHOWN ON ABSCISSA ~ 1 RADIANS/SEC 
MAXIMUM FREQUENCY SHOMKIN ON ABSCISSA 108 RADIANS/SEC 
SrieewmDE LIMITS OF BODE PLOT ARE +-66 DECISELS 
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Figure 9 (Cont.) 


PROBLEM IDENTIFICATION - EX! W’-PLANE T=.01 
XX BODE PLOT CPHASE) XX 


ABSCISSA -> COMMON LOG OF FREQUENCY 
ORDINATE -> PHASE (DEGREES) 

TIC MARKS SHOW MULTIPLES QF 9@ DEGREES 
MINIMUM FREQUENCY SHOLIN ON ABSCISSA = .1 RADIANS/SEC 
MAXIMUM FREQUENCY SHOWN ON ABSCISSA = 1608 RADIANS/SEC 
MAXIMUM PHASE ON ORDINATE SCALE 18@ DEGREES 
MINIMUM PHASE ON ORDINATE SCALE =196 DEGREES 


Figure 9 
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SeeaueENCY RESPONSE 

PROBLEM IDENTIFICATION - EX!1 W’-PLANE T=. 1 

Ee x 0X A AiR KK KK RAKE KAAAR RK RA RAKAAKRKRAKKARAKERELZ 
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Figure 10 w'-Plane Frequency Response Example, T=.1 
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Figure 10 (Cont. )erosLemM IDENTIFICATION - EX1 W/-PLANE T=. 1 


RADIAN FREQ. REAL PART IMAGINARY PART MAGNITUDE PHASE (RAD) PHASE (CEG) 
187226722 ~ 749981617 -44,4216872 46.6277179 -1.58687981 -78,9215483 
» 123284674 ~ 749869949 -48 3467918 48 5537253 -1.58928823 -91.3599392 
141747416 ~. 749828888 -35.2628823 35. 2787736 Siar edo? iz -91,218189 
162975884 ~ 749772759 38 6666494 386758137 ~1.595240863 -91,4085381 
. 187381743 ~ 749699629 -25.6685887 26.6792163 -1.59898859 -91.6182882 
219443469 ~ 749682979 -23.1989122 232030238 -1.68318826 -91.3513782 
1247787636 - 749479252 28 165587 28. 1794298 -1.68794548 S92.\2e0220 
. 284883588 ~ 749386472 - 17 3334525 173494543 -1.61359621 924471289 
3274549 17 749983476 -15.2434888 15.26 18832 -1.61989834 -92 8133714 
376493582 - .748788896 -13.2587881 13.271848 -1.62724549 -93.2343435 
«43287613 - 748399844 -11.5164835 11,3487753 -1.63569823 -93.718128 1 
497782359 - .747886175 -18.8869412 10.8348196 -1,64539465 -94 274283 
1972236769 - 747288246 -3.6925516 $.72446873 -1,465654519 949138817 
697933228 ~ 746313996 ~7.34777687 7 208458359 -1.66935467 955478111 
» 796463332 - 745135214 -§.90026855 §.9925 1428 -1.5848659 34 .4899828 
9469749888 - .743382787 -3..68059274 5 .72985389 -1.78095514 -97 4575856 
1.0880808 | 741548723 492181645 497736475 -1.72833568 -93,5688092 
1.149757 - 738858839 -4.29918654 432279798 ~1.74256117 -99 8414361 
1.32194116 ~ 735343947 -3.6798489 3.75268 159 -1.7689284 -101.398602 
1.5199 189 ~. 738749883 -3.17268863 3.25567824 -1.79717981 ~192.9728889 
1.74752841 ~ 724766284 =eaieer27 dl 2.82237198 - 1, 83849917 -104.879914 
2.88923202 ~.717889 119 S2aGuGlouue 2.444285 14 -1.36851612 -187.859125 
2.31812972 - 707812534 -1.99739681 2.11412199 ~1.91179805 - 189, 53754 
2.588781 ~ .694228885 -1.469848289 1.92555618 - 1.968899 18 =112.301207 
3.85385554 ~ 678841196 ~1.41949479 1.57311962 -2.01641584 Sisisac lay 
3.91119177 -.63779$528 sielaiolo) 1.35211897 -2.97887135 aideliaore 
4,8378173 - 632868475 say idorlore 1,15851483 -2. 14878819 ioe eos 
464158888 ~ 682754245 -. 784821224 . 788939817 -2.22622829 - 127 .553073 
3633669929 ~ 367202992 - .§2833237 340554272 -2.31148516 - 132.4383 1 
§.13598734 ~ 926396326 - .478887285 71181473 248429639 -137.7546885 
7 85489239 ~ 488852854 eacge14ocJ0 397838893 1 -2.08419884 ~ 143.4749 19 
d.1113884 - .431916356 - 293986373 65d 1859494 -2.461880982 - 149 3426 
732683338 SoM LCoe? ~ 178619167 417617384 -2,72072574 139.284 158 
19..7226724 - 330361899 -, 18476951 0346767574 -2.93466462 -162.414377 
12.3284676 ~.28197649 1 ~ 85446999 148 «297233046 -2.94998982 = log eelvoz 
141747418 - 237827592 - 8182671504 237730363 -3,8646771 =lWordvaize 
16.2975886 ~. 196848811 6. 98255476E-83 =. 196960999 3. 18654824 177. 991788 
18.7381745 ~. 161992858 9231913619 » 163644589 2.99939485 171.352796 
21.9443472 -. 132569482 9327930441 13656519 2.8998956 1 166. 186882 
247787439 -, 18827854 8375890315 114616216 2.09749896 158.95744 
28.488359 1 - 8885946878 8398773822 8968299546 2.72618143 156. 198746 
32.7434921 - 872883 1383 9384219724 .8823985323 2.63644429 152.2831 
376493597 ~ 8684938 15 83464598583 978638582 2.999 19126 148.922743 
432876135 ~ 9508148715 8337766115 9618164784 2.59499334 146 .588895 
49 67702345 ~ 9433112758 838770593! 0531298596 2.5238862 144.688879 
27 2236776 ~ 9375272075 8276968923 0466412816 2.00578465 143579936 
$5 .7933236 - 8338884779 8247128768 8412985915 2.50889615 143. 245809 
79 ..6463342 ~ 8296938247 821998753 93698 14454 2.98993127 143.579337 
8.97498 19 - 8271044597 8193382736 83329 13884 2.92207784 144 384467 
100 .880882 prucoicaeveo 8169946678 9383397565 2.9478319 145.93423 
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Peere 10 (Cont. ) 


PROBLEM IDENTIFICATION - EX! W’-PLANE T=. 1 
kX BODE PLOT (AMPLITUDE) Xx 


ABSCISSA -> COMMON LOG OF FREQUENCY 
ORDINATE -> COMMON LOG OF AMPLITUDE 

MINIMUM FREQUENCY SHOVIN ON ABSCISSA .1 RADIANS/SEC 
MAXIMUM FREQUENCY SHOKIN CN ABSCISSA 14@ RADIANS/SEC 
AMPLITUDE LIMITS OF BODE PLOT ARE +-4@ DECIBELS 


Figure 10 
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Peeure 10 (Cont.) 


PROBLEM IDENTIFICATION - EX1 W’-PLANE T=.1 
XX BODE PLOT (PHASE) XX 


ABSCISSA -> COMMON LOG OF FREQUENCY 
ORDINATE -> PHASE (DEGREES) 

TIC MARKS SHOW MULTIPLES OF 90 DEGREES 
MINIMUM FREQUENCY SHOWN ON ABSCISSA = .1 RADIANS/SEC 
MAXIMUM FREQUENCY SHOKN ON ABSCISSA = 108 RADIANS/SEC 
MAXIMUM PHASE ON ORDINATE SCALE = 186 DEGREES 
MINIMUM PHASE ON ORDINATE SCALE = -188@ DEGREES 


Figure 10 
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PROBLEM IDENTIFICATION - EX! S-PLANE 
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Figure 11 s-Plane Root Locus Example 
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Figure 11 (Cont. ) 


PROBLEM IDENTIFICATION - EX1 S-PLANE 
XX ROOT LOCUS PLOT xx 


¢ ¢4444¢44¢4¢4¢44¢4 4 + 


ABSCISSA -> REAL (SIGMA) AXIS 
ORDINATE -> IMAGINARY COMEGAD AXIS 
TIC MARKS SHOW INTERVALS OF i 
THE PLOT FRAME LIMITS ARE: 
OEGSEISSA, -22 TO 11 
ORDINATE, -i? TO 146 
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feeoeenM IDENTIFICATION - EX1 W’-PLANE T=.461 
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ABSCISSA -> REAL (SIGMA) AXIS 
ORDINATE -> IMAGINARY (OMEGA) AXIS 
TIC MARKS SHOW INTERVALS OF 1 
THE PLOT FRAME LIMITS ARE: 
ABSCISSA, -23 TO 12 
ORDINATE, -18 TO 17 


Figure 14 
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C. ROOT LOCUS TEMPLATES IN w'-PLANE 

The s-plane consists of an abscissa representing the 
real portion of the complex variable s and an ordinate 
representing the imaginary part where 

pe Ot JW) 

The w'-plane approaches the s-plane as the sampling 
period T approaches zero. (See [Ref. 11] for a proof.) 
The w'-plane, like the s-plane, consists of an abscissa 
representing the real part and an ordinate representing 
the imaginary part where 

wisu'tj V! 


The s and w' planes are related as shown below. 


Me— 2/1 tanh ( O T/2) 
ee 2/T tan( G)' T/2) 
ee) 2/T tanh" (u'T/2) where -1 < u'T/2 < 1 
Gy = 2/5 tan? Ges 22) where -1 < P'T/2 < 1 


In the s-plane, Q represents damping. Therefore a line 
of constant damping in the s-plane is a vertical line 
parallel to the imaginary axis and passing through the proper 
walue of GO on the real axis. Since the real axis in the 
w'-plane represents u' which is not the damping but only 
meiacted to GC as shown above, it is not possible to easily 
determine the damping of a particular root by Simply ob- 
Serving the position of the root in the plane. This is also 
true of the damped natural frequency, the natural frequency, 


Lae 





ma@aeuhne damping ratio. Because of this it is helpful to 
have templates of constant parameters in the w'-plane to 
make interpretation of the characteristics of a root 
by its location easier. The following sections give some 
insight into the nature of these templates in the w'-plane. 
imeeconstant Damping 

The constant damping templates are shown in Figure 15. 
There are three templates, one of each of three values of 
sampling period, .001 seconds, .01 seconds, and .1 seconds. 
Each template shows lines of constant damping for the values 
25, 590, 75, and 100. These templates are created by using 
the above relationship between O andi u' to find the constant 
meee Of Ww’ that corresponds to the constant value of Q and 
then plotting the constant u' value. From the templates it 
meamseen that for a period of .001 there 1s practically no 
difference from the s-plane. For a period of .0O1 seconds, 
some distortion is noticed as the values of damping approach 
mm@emreciprocal of the value of the period. This distortion 
is in the form of u' becoming smaller than the value of 
damping. And for a period of .1 seconds gross distortion 
is seen. Figure 15 shows an extra template for a period 
oa.) Seconds but for values of damping of 1, 2, 3, and 4. 
In this case only little distortion iS present because the 
values of damping are small compared to the reciprocal of 


eae period. 
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2. Constant Damped Natural Frequency 


Figure 16 shows the templates for constant damped 
natural frequency for values of damped natural frequency of 
25, 50, 75, and 100. The results here are similar to the 
constant damping case. As the period increases, the 
distortion between the value of Gy and the value of aye in- 
creases. In this case the value of VP! becomes greater than 
the value of damped natural frequency as it distorts. 

In the case where the period is .1 seconds the distortion 
is so great that the template is useless and negative 
values of v' are produced by the tangent function in the 
relationship. 

3. Constant Natural Frequency 

Templates for constant natural frequency in the 
w'-plane are shown in Figure 17. In the s-plane constant 
Natural frequency G@), plots as a circle with its center at 


the origin where 
2 2 
Myw= sart( 0 + Gy ) 


To plot constant natural frequency in the w'-plane, substitute 


bey 


ay 


we conn; @ T/2) 


I 


27 ecame Gy) T/ 2) 


into the above equation and solve for Uy ato get 


: 2 a 1 
Vae—eey tan be {Q), - [2/T tanh : carn | 
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For a constant value of natural frequency and period this 
equation is used to plot ut vs Y'. In the w'-plane for a 
period of .001 seconds and a value of natural frequency of 25 
it plots very close to a perfect circle. As the values of 
natural frequency increase they form concentric ellipses 
which are elongated along the imaginary axis and contracted 
along the real axis. This effect becomes more dramatic 

as the period increases. An extra template for a period of 
-O05 seconds is added to better show the trend since the 
change from .01 seconds to .1 Seconds iS So great. 


4. Constant Damping Ratio 


The damping ratio, C , can be defined as 
C= |o Cony | = cos § 


where | represents the absolute value and 6 Tocca e 
formed by the negative real axis and a line joining the 
Origin and the root in a root locus plot. In the s-plane a 
meow Of COnstant damping ratio is a line radially out from 
m@e,O0rigin forming the correct angle - with the negative. 
real axis. To plot constant damping ratios in the w'-plane, 
substitute the proper relations into the above equation to 


get 


aera Cue ye?) 


a, EP: SPC EP I 


{f{tanh > GT /2)4- + [tan_ (p'T/2)]73 


tole 


eS 





Beaving for u' gives 


tan ~ Gay: Tena 


Hee = 2/T tanh 
sqrt(c1i- C* ) 


This equation is used to plot constant damping ratio for 

a given Sampling period. Figure 18 shows templates for 
Semstant damping ratio for values of damping ratio of .1l, 
.o, ./07, and .9. From these templates it can be seen that 
there is a region near the origin that can be interpreted 
the same as the s-plane for a close approximation and that 
the sizeof this region depends on the value of the sampling 
period. Beyond this region a template makes the interpreta- 


tion easier. 
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ABSCISSA -> REAL (UW AXIS 
ORDINATE -> IMAGINARY (NLD AXIS 
TIC MARKS SHCH INTERVALS OF 5S 
WHESPEOGT FRAME LIMITS ARE: 
ABSCISSA, -162 TO S37 
QRDINATE, -168 TO 106 


Figure 15 Constant Damping Templates 
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ABSCISSA -> REAL (UD AXIS 
ORDINATE -> IMAGINARY (NU) AXIS 
TIC MARKS SHOW INTERVALS OF 5 
THE PLOT FRAME LIMITS ARE: 
ABSCISSA, -1464 TO 46 
ORDINATE, -140 TO 198 


Figure 15 
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ABSCISSA -> REAL CLD AXIS 
ORDINATE -> IMAGINARY (NUD AXIS 
TIC MARKS SHOKM INTERVALS OF 53 
THE FECT FRAME LIMITS ARE: 
ABSCISSA, -118 TQ &S2 
ORDINATE, -166 TO 166 


Figure 15 
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Figure 15 
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ABSCISSA -> REAL (U) AXIS 
ORDINATE -> IMAGINARY (NU) AXIS 
TIC MARKS SHOW INTERVALS OF 1 
THE PLOT FRAME LIMITS ARE: 
ABSCISSA, -7 TO 3 
ORDINATE, -5 TO S 
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ABSCISSA -> REAL (U) AXIS 
ORDINATE -> IMAGINARY (NUD AXIS 
TIC MARKS SHOM INTERVALS OF 5 
THE PLOT FRAME LIMITS ARE: 
ABSCISSA, -186 TO 146 
ORDINATE, -37 TO 163 


Figure 16 Constant Damped Natural Frequency Templates 
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ABSCISSA -> REAL (U) AXIS 
ORDINATE -> IMAGINARY (NU) AXIS 
TIC MARKS SHOW INTERVALS OF 5 
THE PLOT FRAME LIMITS ARE: 
ABSCISSA, -166 TO 14¢ 
ORDINATE, -32 TO 168 


Figure 16 
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ABSCISSA -> REAL (UD AXIS 
ORDINATE -> IMAGINARY (NU) AXIS 
TIC MARKS SHOW INTERVALS OF 5 
THE PLOT FRAME LIMITS ARE: 
ABSCISSA, -168 TO 108 
ORDINATE, -195 TG 95 


Figure 16 
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Figure 17 Constant Natural Frequency Templates 
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ABSCISSA -> REAL (U) AXIS 
ORDINATE -> IMAGINARY (NU) AXIS 
TIC MARKS SHOW INTERVALS OF 5 
THE PLOT FRAME LIMITS ARE: 
ABSCISSA, -154 TO 55 
ORDINATE, -165 TO 148 


Figure 17 
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ABSCISSA -> REAL (U) AXIS 
ORDINATE -> IMAGINARY (NU) AXIS 
TIC MARKS SHOW INTERVALS OF 5 
THE PLOT FRAME LIMITS ARE: 
ABSCISSA, -128 TO 77 
ORDINATE, -185 TO 196 


Figure 17 
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ABSCISSA -> REAL (U) AXIS 
ORDINATE -> IMAGINARY (NU) AXIS 
TIC MARKS SHOW INTERVALS OF 5 
THE PLOT FRAME LIMITS ARE: 
ABSCISSA, -126 TO 85 
ORDINATE, -195 TO 100 


Figure 17 
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ABSCISSA -> REAL (UD AXIS 
ORDINATE -> IMAGINARY (NU) AXIS 


TIC MARKS SHOW INTERVALS OF 169 
THE PLOT FRAME LIMITS ARE: 


ABSCISSA, -saaG TO 1408 
ORDINATE, 8 TO 46806 


Figure 18 Constant Damping Ratio Templates 
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Poems 18 (Cont.) 


CONSTANT DAMPING RATIO W’-PLANE T=.81 
XX ROOT LOCUS PLOT xx 
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ABSCISSA -> REAL (U) AXIS 
ORDINATE -> IMAGINARY (NU) AXIS 
TIC MARKS SHOW INTERVALS OF 16 
THE PLOT FRAME LIMITS ARE: 
ABSCISSA, -390 TO 190 
ORDINATE, 8 TO 460 


Figure 18 
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Peeure 18 (Cont. ) 


CONSTANT DAMPING RATIO W’-PLANE T=. 1 
xx ROOT TEGEUS PLOT xx 


S@saeauemeaeae 
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ABSCISSA -> REAL (W AXIS 
ORDINATE -> IMAGINARY (NUD AXIS 
TIC MARKS SHOW INTERVALS OF 1 
THE eeOmeRReME LIM 1 SARE : 
ABSCISSA, —3s4 TO 16 
ORDINATE, @ TO 46 


Figure 18 
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D. COMPENSATION EXAMPLE 

To further demonstrate the use of the RTLOC and FRESP 
programs for analysis and design in the w'-plane, an example 
from [Ref. 6] is duplicated here. 


An angular position servo is described below. 


Md (friction) 
~ + 
GS, ~/) 
LPS ee 
p 





where 
a = 1.9 (V/rad) Ky = 1(A/V) 
3 2 
Ken = 0.385 (Nm/A) Je "Si faces, LO (Nm. Sec) 


Neglecting friction (Md) gives 


Ve 
C, a a 





where 
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Converting to the z domain gives 


9 (2+1) 
G,(z) = k,.T -— 
: : ede 
where 
KK 
K, = —— = 59.8 
2d 


Converting this to the w domain gives 


Ei 1-w 
y= _- 
Z W 
where 
woe= w'T/2 


The open loop transfer function is 
G(w) = Soh aen 
where 


K, = 3.83 


a0 





from the calculation of DC gain. This gives 


9 1-w 
G tw) =—ee, De. > 
W 
For 
© = 0.01 
1-w 
G(w) = 0.0218 > Cay 
W 


Therefore for the uncompensated system equation 5 above is 
entered into the FRESP program producing the output shown 
in Figure 19. From these Bode plots it is seen that the 
uncompensated system is unstable. 
From analog design techniques the final compensation 
network [Ref. 7] is 
One OE we 


K,DCw) = 3.83 oe a 
72 ae) 


and the open loop transfer function becomes 


wtO0.1 
G(w) = K 


w-l 
(w+2.15)° we 
For 


K = -1.017 


ious 





this transfer function was Pees into the FRESP program. 
The resulting Bode plots in Figure 20 show that the compen- 
sated system is stable with a gain margin of 10.5 dB and 
phase margin of 41 degrees. 

Figure 21 shows the output from the RTLOC program for the 
compensated system. This figure includes the normal root 


ieels in addition to an expanded portion of the root locus 


to obtain more detail. These plots are of w and not w' in 
this case. Values of w' can be found by the simple relation 
w' = 2w/T 
where 
T = .Q1 seconds 


for this problem. 

On the expanded root locus plot in Figure 21, a root 
associated with a certain gain may be identified and 
templates similar to the ones discussed in section V.C. may 
be used to easily obtain the characteristics associated with 


the root. 
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Figure 19 Uncompensated System Frequency Response 
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Figure 19 


154 


COSTE TESr Lah = hale 


PHASE (RAD) 


3.83721188 
3.02711786 
3.81687843 
3.90397933 
2.77873683 
2.77839452 
2.968524 

2274338791 
267245986 
prequ4oas 


38177822 


307949031 
33193424 
sodar21 17 
(7302839 
74198151 
60792683 
2.6/844244 
2.63189424 
Dec esr? 

109973266 


PO ro FO A FY ro FS tu 
ses ®@ & a 


2 1428995? 
» 19853559 
d6191371 
$2357 131 
73796531 
70448931 
mcgelocs 
07519926 
036719663 
049218798 
091718948 
(9702562 
of (a2ode2 
76928637 
06438822 1 
(2089989 
2 1437395 
18222842 
1.47846924 
1.698 10376 
1.67946583 


wt em pee OM ee PO OG HP PS OPO PD Pp PO IO PO Oro FO 
= .e = ese 8 = ese &  «@ 


PARMPLE 


PEASE (CEG) 


174.017434 
173.44 1242 
172.593163 
172.13354H2 
elect (a9 
178 929748 
169.525571 
168.957272 
167 254407 
196 .399333 


los. i2uz4 


163.74837? 
162.297198 
149 ..$52798 
133..925748 
197 .d79375 
ido. 18 i200 
193.045 12 
ole ae: 
143.447 103 
145.8435126 
iceocera? 
148. 755479 
133 ..324993 
i159.605278 
Peaecbe(sl 
128 .355593 
127. 749244 
[29,.2802237 
bee eaads 
pea soo70 
113.43/544 
Maa7aloe 
113,.998382 
111.794857 
113.19543! 
ices 40 
106.731905 
W3.c47o02 
184. 228627 
193.0879 15 
101.28672? 
188.9970 16 
D2 aon eos 
79 947197 
Po.cocue co 
gag esoo? 
9§.9692453 
96.2631754 
790¢ 1863 





ieeure 19 (Cont. ) 


PROBLEM IDENTIFICATIGN - KATZ EXAMPLE 
XX BODE PLOT CAMPLITUDE> XxX 


SeeelissA -—> COMMON LOG OF FREQUENCY 
ORDINATE -> COMMON LOG GF AMPLITUDE 
MINIMUM FREQUENCY SHOKIN ON ABSCISSA = .1 RADIANS/SEC 
MAXIMUM FREQUENCY SHOMUN ON ABSCISSA = 16 RADIANS/SEC 
See eUOr LIMITS OF BODE PLOT SRE +-46 OECISEELS 


Figure 19 
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Figure 19 (Cont. ) 


PROBLEM IDENTIFICATION ~ KATZ EXAMPLE 
FA SORE PEOUSCPHASE) xx 


ABSCISSA -> COMMCN LOG OF FREQUENCY 
IROLNATE => PHASE (DEGREES? 

fee Marwce oO MULTIPLES Or 7G DEGREES 
MINIMUM FREQUENICY SHOWN ON ABSCISSA = .1 RADIANS/SEC 
MAXIMUM FREQUENCY SHOKIN ON ABSCISSA = 18 RADIANS/SEC 
MAXIMUM PHASE ON GRDINATE SCALE feeebeGReEs 
MINIMUM PHASE ON ORODINATE SCALE 98 DEGREES 
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Figure 19 
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PROBLEM IDENTIFICATION - KATZ EXAMPLE 
kX BODE PLOT CAMPLITUDE) Xx 
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ABSCISSA -> COMMON LOG OF FREQUENCY 
ORDINATE -> COMMON LOG OF AMPLITUDE 
MINIMUM FREQUENCY SHObIN ON ABSCISSA = .1 RADIANS/SEC 
MAXIMUM FREQUENCY SHOBIN ON ABSCISSA = 14 RADIANS/SEC 


Ciel yUOS LIMITS OF BODE PLOT 


Apevet= od DEGCTEELS 


Figure 20 
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PROBLEM IDENTIFICATION - KATZ EXAMPLE 
AA SEODE PLOT <PHASE. kz 


ABSCISSA -> COMMON LOG OF FREQUENCY 
ORDINATE -> PHASE (DEGREES) 

Hic MARKS SHOL! MULTIPLESSOF 4a DEGREES 
MINIMUM FREQUENCY SHOWN ON ABSCISSA = .1 RABDIANS/SEC 
MAXIMUM FREQUENCY SHOVIN ON ABSCISSA = 10 RADIANS/SEC 
MAXIMUM PHASE ON ORDINATE SCALE = 186 DEGREES 
MINIMUM PHASE ON ORDINATE SCALE = -18@ DEGREES 


Figure 20 
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Figure 21 Compensated System Root Locus 
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PROBLEM IDENTIFICATION - KAT2 EXAMPLE 
Ae ROS TEREUS FEO) xx 


ABSCISSA -> REAL (SIGMA AXIS 
ORDINATE -> IMAGINARY (OMEGA) AXKIS 
TIC MARKS SHOW INTERVALS OF 1 
THE PEO) eRe Ei ite eRe: 
ABSCISSA, -8 TO 2 
ORDINATE, -S TO 3S 


Figure 21 
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PROBLEM IDENTIFICATION - KATZ EXAMPLE 
ye ROOT LOCUS PEOT Ax 





ABSCISSA -> REAL (SIGMA) AXIS 
ORDINATE -> IMAGINARY (OMEGA) AXIS 
TIC MARKS SHOW INTERVALS OF 1 
THE PLOT FRAME LIMITS ARE: 
ABSCISSA, -1 TO 2 
ORDINATE, -1!1 TO 2 


Figure 21 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


This thesis demonstrates that common classical control 
programs may be adapted to run on an inexpenSive microcomputer 
system. With the exception of the memory management problem 
encountered, this was done with relative ease. This allows 
more people access to these programs by introducing a new 
group of computers on which they can be run. Also these 
programs require no knowledge of any computer language or 
input card formats making them easier to use than the 
Fortran versions. 

It is also demonstrated that the transfer function pro- 
grams are a useful tool in the study of sampled data systems 
as well as classical systems. 

Only a sampling of five programs was converted in this 
thesis. Other programs can and should be Similarly converted 
to run on microcomputer systems. There are programs in 
existence for the same microcomputer system used in this 
thesis that generate the aircraft stability derivatives. 

An effort should be made to modify these programs and the 
programs of this thesis to make them compatible and 


complementary. 
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APPENDIX A 


DESCRIPTION OF MICROCOMPUTER SYSTEM 


The microcomputer system used in developing this thesis 
consisted of the following components: 

Apple II plus computer (48K) 

Sask It 3 1/4" floppy disk drive and controller card 

USI 9" green screen monitor 

NEC PC-8023A-C dot matrix printer 

Grappler printer interface card 

Add Ram 16K expansion card 

The programming language used was Applesoft basic. All 
graphs were generated using the High Resolution graphics 
commands. In the High Resolution graphics mode a 
matrix of dots 280 dots wide and 192 dots high can be dis- 
played. The High Resolution page one is the only page used 
in this thesis and it resides in memory in the 8,192-byte area 
from $2000 to $3FFF. High Resolution page two resides in 
the area from $4000 to $6000. 

A memory management problem was encountered during the 
programming of the programs requiring the use of the graphics 
capabilities of the microcomputer system. To understand the 
nature of the problem encountered a brief description of the 
normal use of memory in the Apple II is necessary. An Apple- 
soft program is normally loaded at memory location $800 and 
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loads up. LOMEM is set to the end of the program. Simple 
variables are stored from LOMEM upward as they are defined 

in the program. This gives only 6K-bytes of RAM before there 
is a conflict with the first High-Resolution page and only 14K 
before there is a conflict with the second High-Resolution 
page. [Ref. 8] Since the programs and variables in this 
thesis exceed 14K, simply using High-Resolution page two in- 
stead of one is not the answer. The first part of the solu- 
tion was to get the Applesoft programs to load above the 
Space used by High-Resolution page one and in effect pro- 
tecting that space from interference. This fix creates a 

new problem. The disk operating system loads at the top of 
the 48K of memory and sets HIMEM to $9CF8. String variables 
start at HIMEM and build down. This allows 22K of space for 
the program and all variables. Some programs in this thesis 
require more space than this. Even with a 16K memory 
expansion card installed the disk operating system will 
ignore it so the added memory is useless. This final problem 
was solved by uSing a utility program that relocated the disk 
operating system into the higher memory provided by the 
memory expansion card and resets HIMEM to $BFOO. This gives 
32.9K of useable program space which is sufficient for all 


programs in this thesis. [Ref. 9] 
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